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The goal of this paper is to exhibit a deep relation between the partition function of the Ising 
model on a planar trivalent graph and the generating series of the spin network evaluations on the 
same graph. We provide respectively a fermionic and a bosonic Gaussian integral formulation for 
each of these functions and we show that they are the inverse of each other (up to some explicit 
constants) by exhibiting a supersymmetry relating the two formulations. 

We investigate three aspects and applications of this duality. First, we propose higher order 
supersymmetric theories which couple the geometry of the spin networks to the Ising model and for 
which supersymmetric localization still holds. Secondly, after interpreting the generating function 
of spin network evaluations as the projection of a coherent state of loop quantum gravity onto the 
flat connection state, we find the probability distribution induced by that coherent state on the 
edge spins and study its stationary phase approximation. It is found that the stationary points 
correspond to the critical values of the couplings of the 2D Ising model, at least for isoradial graphs. 
Third, we analyze the mapping of the correlations of the Ising model to spin network observables, 
and describe the phase transition on those observables on the hexagonal lattice. 

This opens the door to many new possibilities, especially for the study of the coarse-graining and 
continuum limit of spin networks in the context of quantum gravity. 
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I. INTRODUCTION 


This paper is devoted to studying the relationship between two fundamental objects in physics, associated to finite 
graphs: the two-dimensional Ising model and the spin network evaluations on planar graphs. 

Given a graph T and a coloring of the edges of F, i.e. a map c : E —► N where E is the set of edges, the 
classical spin network evaluation (F, c) is a rational number which is the result of contracting some tensors over 
irreducible representations of SU(2). Spin networks arise in many areas, in particular related to physics. Since they 
come from the representation theory of SU(2), they are objects of prime interest in the theory of quantum angular 
momentum [li where they are often called Wigner symbols. As such, they have applications in atomic/molecular 
physics, chemistry, quantum information and so on J2J . More recent applications stem from quantum gravity as 
spin network equipped with holonomies are the states of loop quantum gravity while their evaluations provide 
quantum gravity amplitudes, known as spin foams, m- The latter are intimately related to lattice topological 
invariants, such as the Reidemeister torsion (ehei and the Turaev-Viro invariant of 3-manifolds. 

Spin network evaluations can be computed in many different ways. The most famous is certainly the combinatorial 
definition due to R. Penrose mi- In quantum gravity one often uses contractions of SU(2) intertwiners (notice that 
it requires an orientation on the edges while Penrose’s definition does not - we will prove the equivalence between 
those evaluations in the main text). In the last years it has been understood [1311131 that a good way to study spin 
network evaluations on a fixed graph T is to organize it in a single generating series Z spin (T) := J[] c (r, c)Y c where 
the symbols Y c stand for a suitable multivariate monomial in formal variables Y e ,e € E (full details will be provided 
later). 

On a seemingly different side of physics (and mathematical physics), the 2D Ising model can be defined on the 
same graph F. It probably is the most famous statistical model, based on the configuration space of maps from the 
set V of vertices of T into {±1}- The Hamiltonian (energy function) of the model is a sum of interactions between 
nearest-neighboring sites of T, hence associated to the set E of edges of T, and weighted by couplings y e , e £ E. The 
partition function Z Ising (T) was proved by van den Waerden to be proportional to a sum over even subgraphs of T 
weighted by some monomials in tanh(j/ e ),e £ E(T) (see [T7] for instance). 

The present paper is motivated by the following observation made in EU; if r is a planar trivalent graph and for 
each edge e £ E we set Y e := tanh(y e ) then the following equality holds: 

(^Z Isin 3 = ^n cosh (ye)j 2# 2V . (1) 


Such a relation was independently noticed and shown to hold for the square 2D lattice with homogeneous couplings 
in 122) . A straightforward proof of that equality can be given using van der Waerden high temperature expansion for 
the Ising model and Westbury’s formula for the generating series of spin networks. However, this approach sheds no 
light on the intimate reason why the equality holds. 

In the present paper we explain this phenomenon by identifying a supersymmetry which relates the Ising model on 
a trivalent planar graph T to the generating series of spin networks on it. In order to achieve this, we: 


1. represent Z Ismg and its square via fermionic Gaussian integrals on a space of fermions indexed by the half-edges 
of r (Propositions II.3 and II.5), 


2. compute Z spm as a standard (bosonic) Gaussian integral on C 1 2 3 # E (Theorem III.5), 

3. provide a supersymmetry relating the fermions and bosons associated to each half-edge (Section |IV A). 


Formulas for Z Ismg based on fermionic (or Grassmanian) integrals have existed since the beginning of the 80’s (see 
[23], [24] , [25j 1. Our contribution here is to provide a new formula in which the integrand is a Gaussian whose bilinear 
form is defined using a Kasteleyn orientation on F (which always exists by [26J, [27! as F has an even number of 
vertices). Our formula is inspired by the integral representation of [23! on the hexagonal lattice (where Kasteleyn 
orientations are not required though). 
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On the spin network side, a first Gaussian integral representation for Z spm was found in m (see also dH |TS] for 
another Gaussian integral approach applicable for non-trivalent graphs), where a regularization process was needed 
to compute the integral. Here we provide another (very similar) formula via a convergent integral but based on a 
space whose dimension is twice as large. Another observation here is that the definition of Z spm we use depends on 
an orientation on F while standard spin networks do not. In Theorem III.4 we prove that our evaluation still coincides 
with the standard one, provided the orientation is a Kasteleyn orientation. 

The main point in the above formulas is that they really “look similar”. This is formalized in Section |IV A| by 
means of a supersymmetry between the fermionic and bosonic degrees of freedom. We show that a Berezin integral 
whose argument is the product of the above two integrands equals a constant, by a simple supersymmetric localization 
argument. 

This unveils a deep relation between two precedently unrelated objects. We start exploring aspects and consequences 
of this in Subsection |IVB and then in Section [V] In Subsection IV B| we introduce a supersymmetry-preserving 
generalization (whose extensive study will be pursued elsewhere) where the Ising model and the spin networks on T are 
coupled non-trivially. In particular, the Berezin integral is not Gaussian anymore. We argue that the supersymmetric 
localization we have found will be useful to compute such modified versions of the generating function of the spin 
networks evaluations. 

In Section [vj we start by recalling the origin of the generating series Z spm in loop quantum gravity, as a coherent 
state (see also [Hj and m) whose background geometry is set by the couplings Y e and where the spins on the 
edges are quantum numbers of length. The expectations of products of length operators in this coherent state are 
moments of a probability distribution on the spins which depends on the couplings Y e , e € E. The stationary phase 
approximation of this distribution, at large spins, leads to relations between the couplings Y e and the values of the 
spins at the stationary points. Using the physical meaning of the spins as quantum numbers of length, those relations 
can be written in geometric terms. Remarkably, they are found to be the same as those which define the critical 
couplings of the 2D Ising model on isoradial graphs E3E0]. In particular, the spins only determine the geometrical 
shape at the stationary points and can be arbitrary rescaled. 

the fundamental equality 0 between Z spm and Z Ismg to relate the observables 


We further use, in Section 


VB 


of the Ising model, i.e. Ising spin correlations, to length observables in loop quantum gravity, in particular to matrix 
elements of products of length operators between the coherent state and the physical state (i.e. the flat connection 
state). In Section VC we push this further to compute the generating series of the matrix elements of the length 
operator of a single edge to an arbitrary power: it takes a simple closed form as a function of the nearest-neighbor 
correlation function of the Ising model. 

Finally in Section m another integral representation for coherent spin network states is introduced. It is based on 
integrals over spinors (variables on C 2 ) associated to the vertices of the graph. This way, we show that one can define 
other generating functions which differ from Z spm by choices of combinatorial factors whilst still admitting Gaussian 
integral representations. We find in particular an example of such generating functions for which the stationary points 
of the distribution on the spins for the matrix elements of length operators between that coherent state and the flat 
connection state depend on the scale of the geometry, in contrast to the result obtained for the ordinary generating 
function Z spm . This generalizes the comparison between different choices of generating functions of spin network 
evaluations started in [14] on 2-vertex graphs. 

The questions opened by the present work are multiple and we hope that the supersymmetry we here introduce 
between the 2D Ising model and the spin network evaluations will allow to exchange and cross-fertilize results in 


Statistical Mechanics and Loop Quantum Gravity. We outline some of those perspectives in Section VII 


Acknowledgements The three authors profited of a PEPS funding from CNRS for the project “Spin-Ising”. 
Notations. In all the paper we will let T be a finite graph, V be the set of its vertices, E the set of its edges, E the set 
of its half-edges or, equivalently of its oriented edges (E has a natural 2 —>■ 1 map to E), and A the set of its angles i.e. 
pairs of distinct half edges having the same endpoint. For each oriented edge e : s(e) —>• t(e), we write s(e),t(e) £ V 
for its source and target vertices. 

Unless explicitly stated the contrary we will assume that F is planar, i.e. embedded (up to isotopy) in S 2 . This 
automatically equips T with the datum of a cyclic, counter-clockwise ordering of the edges around each vertex. For 
each angle a around a vertex, following the cyclic ordering around that vertex, we call s(a) the source half-edge of 
the angle and t(a ) its target half-edge. 

Reciprocally, a cyclic ordering of the edges around each vertex allows to canonically thicken F to an oriented surface 
with boundary which we will denote S(r): since F is assumed to be planar, S(T) is homeomorphic to the complement 
of a collection open discs in S' 2 . A connected component of the complement of S(F) is called a face of F and it is 
naturally equipped with the counter-clockwise orientation. 

F will further be assumed to be connected and bridgeless (i.e. 1-particle irreducible, or, equivalently no edge 
disconnects T). Taking bridges into account is quite simple but requires to extend some definitions (like the equivalence 














4 


class of Kasteleyn orientations) and it does not bring much to the theory. Recall indeed that from recoupling theory 
the spin of a bridge in a spin network has to vanish, which thus factorizes the spin network evaluation into two parts 
associated to disjoint graphs. Note that being bridgeless is equivalent to saying that each edge is incident to exactly 
two faces. 

In the following we will associate functions to T which depend on one of two kinds of parameters: the Y e variables, 
indexed by e £ E and the X a variables, indexed by a G A. 


II. 2D ISING MODEL 
A. Loop Expansion of the Ising Model 

Definition II.1 (Ising Model). An Ising spin configuration on T is a map cr : V —> {—1,+1} associating ±1 to each 
vertex of the graph. The partition function of the Ising model on T is a function of couplings y e along each edge: 

Z Ism9 (T, {y e }) = exp | ^2 yecrs(e)v t (e) 
cr VeGE 

The van der Waerden identity exp(ya u a v ) = cosh(K)(l + tanh(y)a u a v ), for all a UtV = ±1, allows to re-express the 
partition function as 

Z Islng (T , {y e }) = ( U cosh(y e )) ^ J|(l + tanh(y e )cr s(e) er t(e) ) 

eGE a e 

= 2 # v (n cosh(y e )) E ntanh(y e ) = 2 #V (]J cosh(y e )) En^ 

eGE eG7 eGE 7(E£/ eE7 

where we write Y e = tanh(y e ) and we sum over the set Q of even subgraphs of T (also known as Eulerian subgraphs) 
i.e. subgraphs 7 C T such that every vertex of G is incident to an even number of edges of 7 . 

In this paper, we will focus on 3-valent graphs, i.e such that every vertex has exactly 3 edges attached to it. If T is 
a 3-valent graph, then all even subgraphs are unions of disjoint loops. This is the high-temperature loop expansion 
of the Ising model, which we will match against the evaluation of spin networks. 



B. Grassmannians for the 2D Planar Ising Model 

Definition II.2 (Kasteleyn orientation). A Kasteleyn orientation on T is an orientation of the edges such that each 
face has an odd number of edges whose orientations do not match the one induced by the face. 

When drawing T on the plane, it means each face has an odd number of clockwise edges, see Fig. [T] If T is 
embedded in the plane (and therefore the outerface, i.e. the connected component of its complement that is infinite, 
is not considered as a face), then Kasteleyn orientations exist (Kasteleyn’s theorem). If T is embedded in S 2 , then the 
notion of outerface is meaningless as which face is the outerface in the drawing depends on the choice of projection on 
the plane. If one insists on drawing T in the plane, the outerface then has to have an odd number of counter-clockwise 
edges. It is known that Kasteleyn orientations on cellular decompositions of oriented compact surfaces exist if and 
and only if the number of vertices is even EH, EZj. A regular graph of degree 3 has an even number of vertices. 

We introduce a set of Grassmann variables attached to half-edges, (V’s(e); V’i(e) }e. which all anti-commute with one 
another. We define an edge action for each edge e and a corner action for each angle a: 


Ie(lPs(e)i' l Pt(e)') tA(e) ^t(e) > At ilfs(a) > 'At(ct)) I^s(ct) t^i(ct) ■ (2) 

In what follows we shall also use the notation iff, to denote the Grassmann variable associated to the half-edge 
contained in e and incident to v (if it exists it is unique as T is bridgeless). 

Proposition II.3 (Grassmannian expression of the Ising model). The partition function of the Ising model on a 
3-valent planar T (equipped with a Kasteleyn orientation) reads 

Z Isin 9(r,{y e })= (jJcosh(y e )j 2 # v Z f (T,{X a }) , 


( 3 ) 
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FIG. 1. Example of a Kasteleyn orientation around and inside a cycle c (in bold): each face of the planar graph contains a 
odd number of clockwise oriented edges on its boundary, which implies that the number of clockwise oriented edges around an 
arbitrary cycle is related to the number of interior vertices; here F = 3, V m t = 1, Ei n t = 3, a(c) = 3 (a is a large angle). 


Z f (T,{X a }) = f n #t(e)# S (e) exp ( X! /e ('^(e)’^(e)) + Y1 X » 


where X a = (y 5 ( a )lf( a )) 1 ' /2 . The Grassmannian integral is normalized as f Y[ e dift(e)d4’s(e) Tie '0*(e)'0»(e) = 1- 

The proof is based on the following lemma on Kasteleyn orientations. Let c be a cycle of T. When T and c are 
drawn in M 2 , there is a well-defined inside and outside of the cycle. The disc inside c induces a counter-clockwise 
orientation on c. A large angle of c is a pair of half-edges of the cycle both incident to the same vertex and such that 
the third half-edge incident to this vertex (and not contained in c) lies on the inside, as illustrated in fig(T] 

Lemma II.4. Let c be a cycle, a(c ) its number of large angles, #E c ;(c) its number of clockwise edges, and ff\li n t{c) 
the number of vertices on the inside. Then 

(-l)“(c) = (_ 1 )#v i nt(c) j ( 4 ) 

(_!)#Ecl(c) = + (5) 

Proof. We begin with the equation Q. Let F(c) be the number of faces composing the disc bounded by c and #E int 
be the number of internal edges in such disc. The subgraph made of the region inside the cycle c and bounded by 
it (erasing all vertices in c where c does not form a large angle) is planar. Computing the Euler characteristic of S 2 
cellularized via this planar subgraph we get: 


F(c) + 1 + (i(fV int (c) + a(c)) — (#E int{c) + a(c)) — 2 . 


( 6 ) 


The extra +1 counts the external face to the cycle c. Moreover, the number of vertices on the boundary cycle c 
obviously equals the number of boundary edges, thus we are left with counting only internal vertices and internal 
edges. On the other hand, we translate the fact that the graph is 3-valent into an equation relating the number of 
internal vertices and edges: 


3#(V int (c) + a(c)) = 2(#E int (c) + a(c)). (7) 

Combining these two equations yields: 

a(c ) = 2(F(c) - 1) - #V mt (c), (8) 

which implies the desired result (El). 

We now prove the equation ([S]). The quantity on the left hand side is obtained by assigning the weight —1 to 
clockwise edges and 1 to counter-clockwise edges of c and taking the product of those weights. We describe another 
way to do that. One proceeds the same way on each of the F(c) faces in the interior, to get weights e/. The product 
of those weights is not exactly what we want: each edge comes once with a +1 and once with a —1. Thus, it has to 
be corrected via a factor (—l)# E i"t( c ). Thus we get 

(_l)#Ed(c) = JJ e/ x 
/ 


( 9 ) 
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Moreover, since the Kasteleyn orientation ensures an odd number of clockwise edges to each face, e/ = (— 1) F ( C ). 
Using again Euler’s relation <©> one gets 


(_l)#Ed(c) = (_i)F(c) (_ 1 )F(c)-l-#V int (c) j 


( 10 ) 


II.4 


which simplifies to ([5]). 

Proof of the Proposition \ II. 4 We just have to show that the fermionic integral equals the expansion onto even 
subgraphs of T, 


n #i(e)#s(e) ex P I X a^s{a)^t{a) ) = II Ye 

e \ e a / 7 £C? eG 7 


(ii) 


for X a = (y s ( a ) y^Q,)) 1 / 2 . Notice that 7 € £/ being even, and T being regular of degree 3, 7 is a disjoint union of 
cycles. 

We simply expand the exponential of each edge and corner actions within the integral and commute the sums with 
the integral. Each exponential terminates at the linear order and only terms which saturate the number of ips survive. 

If a term comes from I a for some corner a made of the half-edges h , g, then I e cannot appear in this contribution 
if e is an edge containing h or g. Therefore the Grassmannian variable incident to g (or h) and contained in the same 
edge as g (or h) must come from a corner action and so on, until it closes to a cycle c. A typical contribution is thus 
labeled by a collection of disjoint cycles, i.e. an even subgraph 7 £ G, as expected. All half-edges not contained in 7 
must come with their corresponding half-edges to form full edges and they get from I e the weight 1. Moreover each 
cycle receives the weight n a 6 c*a = n e6 c Y e , up to a sign which is crucial to determine. 

Consider a cycle c and let #E(c), #E c ;(c) be respectively the number of edges c is formed of and the number of 
edges in c which are oriented in the clockwise direction (here use use the fact that c bounds a unique disc in R 2 ). 
Since all the terms I a commute with each other, we can reorder their product freely. We label the vertices of c from 1 
to I following the counter-clockwise orientation and order the I a accordingly; let us denote by ip* the Grassmannian 
iph where i denotes the vertex the half-edge h is incident to and j the vertex such that h belongs to e : i -H- j. The 
product of Grassmannians we get from expanding the exponentials is then 


(-1 )#B(c)-a(c) ^1^2 . . . = (•_ 1 )#E(c)-a(c) + l 


tp{ip\ipltplipl ■ ■ ■ ip e e_ 1 


( 12 ) 


where we reordered the variables ip s ( a )ipt(a) exactly for the small angles i.e. where their order did not match that 
induced by the orientation of c and we have brought ip{ to the beginning so that the Grassmannians which lie on the 
same edge are next to each other, in the counter-clockwise order. Next, in order to match the measure Y[ e ^V’t(e)^s(e) 
we flip the order of those Grassmannian variables associated to half-edges which belong to an edge oriented counter¬ 
clockwise. The remaining integral evaluates to 1. Therefore the sign of the cycle c is 


(_l)# E ( c )- a (c) + l x c! (c) _ 


(13) 


by Lemma II.4 As a conclusion, all cycles have a positive sign. 


II.4 


We can further push this path integral reformulation of the Ising model by using complex fermions instead of real 
fermions. This is a non-essential mathematical trick, doubling the number of variables associated to each half-edge of 
the graph (thus having other anticommuting variables iph), but it will later allow to write explicitly the supersymmetry 
relating the Ising model to the spin network evaluations: 


Proposition II.5 (Complex Grassmannian expression of the Ising model). The partition function of the Ising model 
on a 3-valent planar T (equipped with a Kasteleyn orientation) reads 


z Ising (r, {ye}) = Y[ cosh (y e ) 2 # V Zf(T, {AJ) . 


(14) 


VeeE 


zf(r,{x a }) = J ll #t(e)#s(e)#i(e)#s(e) exp 4>s(e)A(e) 1ps{a)^t{a) j , 


where the Grassmannian integral is normalized as / lie dipt(e)dip s (e)dipt(e)dip s (e) Tie l/ ; t(e)V’s(e)'^t(e)' ! /’s(e) = 1; and we 
denote ip l ( or ipf a variable associated to the half-edge contained in an edge e and incident to the vertex v. 
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Furthermore it holds: 

h 


Zf(T,{X a }) 2 = J 


This can be proven by either repeating the same steps as above with the detail of orientations and signs, or by 
directly performing the Grassmannian integrals over the variables ip”. We detail below the explicit computation of 
the integral tracking all the signs. 


Proof. The proof of the formula for Z9 is similar to that of Proposition 


II.3 


we need to find all the ways to obtain 


the top monomials YlevW by taking products of monomials of the forms ip”ip”, ( — Vye)V’t(e)) and ipstaVua) and 
for each such way we need to check that the overall sign of the integral is 1. 

Observe that in each such way of factorizing the top monomial, each time a factor ip s r a )ipt(a) is present then 
necessarily also a factor (—Vb(e)V ; i(e)) is present with e being the edge containing t(a). Similarly, each time a factor 
(■ —ipsieVtte )) is present then so is a factor of the form ip 8 ( a )ipt(a) where a is one of the two angles at the end of e. 
This implies that each factorization of the top monomial corresponds to a disjoint union of cycles c (corresponding to 
closed loops of factors Yl aec if s ( a )4’t( a ) TLec VMe)^t(e)) and a monomial of the form IL^c i>*(e) , <p8(e)ipt(e)‘>Pt(e)- 

Now orient each component of c as induced by the disc it bounds in R 2 and if ex,... e*, are the edges encountered 
while circulating in c and Vi = e^_i H ej (iq = e*, fl e 4 ) reorder the terms ip” and ip” to get 


±Ve\reirel^ 


■VetreXr 1 


e^c 


e)Vb(e)^t(e)^(e) = ±{Ve\Vel ' ' ' VetVe h ) (' Ve \# 


Kl 


■VelV 1 


Zk )n Vb(e) Vb(e) V*t(e) ^Pt(e) ■ 

e^c 


Let us analyze the sign we get from the right hand side. We permute the variable iph exactly as in the proof of 
Proposition II.3 we first reorder internally all the degree two monomials ip s ( a )' l Pt(a) according to the order induced 
by c, thereby acquiring a —1 sign if and only if a is a small angle. Then we permute ipj. with all the other terms ip], 
thus acquiring a —1 sign. At this stage we get the monomial: 


(_l)#E(c) a(c)+i(^1^2 . . . ^l)(H -1p s {e)1pt(e))([[ VV) = (~l) #E(c) a(c)+1 (VqVl • ‘ • ^l)([{Me)V{e)){\{ W) 

e€c h£c eSc h^c 

= (-l) 2 # E( c ) -# E c , + 0 ( c) + 1 (n^(^ s ( e ))( n ^(^ s(e ))(n^^) = (Y[A(e)V(e)V(eV S (e)){Y{Mh) 

h£c eGc h£c 


eGc 


eGc 


where in the second equality we reordered the terms iplipl +1 if and only if the edge was oriented counter-clockwise 
(i.e. if and only if —> Vi + 1 ), and in the third we used Lemma II.4 So the Grassmannian integral of this term is 

+n aC c*« = Ileecbe as claimed. 

Before moving to the proof of the second statement, let us first define antisymmetric matrices A and B of size 
2 #E x 2 #E as follows: 


Y -^s(e)^(e) = \V 


y ^X a ip s(a) ip tM = Bip. 


Then the first statement implies the following: 


Zf(T,{X a }) =±Pf 


B 

Id 

-Id 

A 


±\/det (I + AB) 


where Pf(X) = Ak /det(X) is the Pfaffian of the antisymmetric matrix X , in the last passage we used the fact that the 
determinant of a 2 x 2-block matrix (mImI) det(M 4 — M 2 Mf 1 M 3 ) det(M 4 ) if M 4 is invertible and here M 4 = A 
is such that A 2 = —Id. 

We now remark that the integral considered in the second statement equals the Pfaffian of the antisymmetric matrix 
M of size 8#E x 8 #E which in the basis given by {ip e }, {he], {V’e}, {h e } i s: 


( B 

0 

0 

-Id \ 

0 

B 

-Id 

0 

0 

Id 

A 

0 

V Id 

0 

0 

A J 


M = 


(16) 

















Observing that A 2 = — Id and using again the fact that the determinant of a 2 x 2-block matrix ( ^ ^ ) is det(Mi — 
M 2 M 4 - 1 M 3 )det(M 4 ) if M 4 is invertible, we get: 

det(M) = det(A) 2 det(A + B) 2 =$■ Pf(M)=det(A + B) = Zf{T,{X a }) 2 . (17) 


II.5 


C. From Edge Variables to Angle Variables and Back 


Although we started with the Ising partition function Z Ism9 (T,{y e }) as a function of the edge variables y e (or 
Y e = tanhy e ), it was convenient to switch to angle variables X a to define and compute its reformulation as an 
odd-Grassmannian integral. This relies on the following fact: 

Lemma II. 6 . Given variables Y e on the edges of a graph T, assumed to be planar, connected and bridgeless, if we 
define angle variables X a = (Ts(a)^t(a)) 5 > then we have the following equality for any closed loop 

n =n y e . (is) 

ae£ ee£ 


It is interesting to invert this mapping and check if the equality between the Ising partition function and the 
fermionic path integral holds for arbitrary angle couplings X a . Of course, the number of angle variables on a 3-valent 
graph is twice the number of edge variables. It is nevertheless possible to reverse the above mapping: 

Lemma II. 7. Given variables X a on the angles of a graph T, assumed to be 3-valent, planar, connected and bridgeless, 
if we define edge variables following the convention of the left part of Figure^ 7} 

... _ f X GG1 X ee2 X eGl X eG2 

e l Y Y- - 

\ y '-e 1 e 2 y '-e 1 e 2 



where e 4:2 are the two other edges incident to the source vertex s(e) while e 4i2 are the two other edges incident to the 
target vertex t(e), then we have the following equality that holds for any closed loop C: 

n >■ = n 

ee£ a££ 


These two mappings allow to 
We further notice that if we 
mapping (191, then we recover 
mappings are not one-to-one. 


consider as primary variables either the edge variables Y e or the angle variables X a . 
start with edge variables Y e and we apply the relation X a = (T s (a)^(a)) 5 to the 
the initial variables Y e . Clearly the reverse is not true, due to the fact that these 


III. THE GENERATING FUNCTION OF SPIN NETWORK EVALUATIONS 

We now turn to the realm of quantum geometry and 3d (Euclidean) quantum gravity, where the quantum states 
of 2D geometry are constructed from the 3nj-symbols of spin recoupling PUSH, and so are transition amplitudes for 
instance in the Ponzano-Regge model (see 7], and mmm for more recent developments). 


A. The Hilbert Space of Spin Networks and their Evaluations 

We consider a planar, 3-valent and connected graph T. We equip the graph with a counter-clockwise cyclic ordering 
of the edges incident to each vertex. Let us also choose an orientation o of the edges of T. For the correspondence with 
the Ising model (which requires to get a sign + for all loops), this orientation will be chosen to be of the Kasteleyn 
type. 

We will discuss the (in) dependence of the spin network evaluation with respect to the chosen orientation and the 
relation between our (tensorial) definition of evaluation and the other standard definitions of 3nj-symbols in the next 
section. 
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In the context of loop quantum gravity and spinfoam models P3H> we define the Hilbert space of wave-functions 
for the 2D quantum geometry on the graph T as the space of gauge-invariant functions of SU(2) group elements living 
along the oriented edges of the graph: 

/({SeW) = fdKfogehsuher), V{/i„}„ eV G SU(2) #V , (20) 

whose scalar products are defined by the Haar measure on SU(2): 

= f [dg e \ 4>({ge}) 4>{{9e\) ■ (21) 

JSU( 2 )# E 

A basis of this Hilbert space L 2 (SU(2)# E /SU(2)# V ) of square integrable functions over SU(2)# E is provided by the 
spin network functions. For an oriented 3-valent graph T, spin network states are labeled by a graph coloring. 

Definition III.l. A coloring of T is a map col : E — > N; we will denote the value of the coloring on an edge e by 2j e , 
and we will call j e the spin of the edge for the given coloring. 

The spin characterizes an irreducible representation of SU(2) of dimension (2 j e + 1) whose basis vectors are indexed 
as e J m (also denoted \j,m)), for m £ {— j e , —j e + 1,..., j e }; these vectors form the magnetic number basis. The spin 
network function is then defined as: 


Eel 


(».) = e Ih-v 

{ ml} e 


D :h t 

—ml 


Me) n 


Je% Je 


mz, 


mz 


m'z 


( 22 ) 


where the edges around the vertex v are cyclically ordered as j e «, , j e ;> with the magnetic moment labels m’ living 

on every half-edge around every vertex. Finally ) = (j, m\g\j, m') is the Wigner matrix representing the 

SU(2) group element g in the irreducible representation of spin j. The symbol (yA h h 'j j s the Wigner 3j-symbol 
(see [I], Chapter 3.7). It is invariant under rotations, meaning that it is an intertwiner between the tensor product 
of the three representations and the trivial one. In other words, the vector J 2 m i m 2 m 3 (m) m 2 m 3 ) e m i ® e m 2 ® e m 3 
is a basis of the (1-dimensional) space of SU(2)-invariant vectors in the tensor product of the three representations. 
The entries of the projector onto the invariant space can be written as 



dgD^ ni {g)D^{g)D%X 3 {g) 


Jl J2 


, 777-1 7772 TO 3 



Jl J2 J3 
n\ n 2 n 3i 


The orthogonality and completeness relations read 


E 

mi ,771-2 


Jl J 2 J 
mi TO 2 to - 


[ Jl J 2 f 

l TOi TO -2 TO' 


1 

2 j + 1 


djj' 


S 


mm' 1 


E(^' +1 ) 

j,m 


Jl J 2 J 

mi TO 2 to. 


f Jl J 2 j 
l rn\ m' 2 m 


d m2rn ^. 


Wigner 3j-symbols are invariant under cyclic permutations of their three entries and pick up a sign upon reversing 
the cyclic orientation as well as upon reversing the orientations of the incident edges: 


Ji J 2 J 3 

- TO -1 TO 2 TO -3 


= (- 1 ) 


Jl+j2+j3 | Jl J3 J 2 
\TOi TO 3 TO 2 : 


= (-l) 


71+72+73 [ Jl J 2 J3 

\-TOi TO -2 -TO 3 


= < 5 , 


'mi+m2+m 3 ,0 


Jl J 2 J3 

,TOi TO 2 TO 3 


The central object of our paper is the spin network evaluation, defined as the evaluation of the spin network function 
at the identity on all edges, g e = I. 

Definition III. 2. For a colored 3-valent graph T, equipped with an edge orientation o, we define its “unitary (up to 
sign) spin network evaluation” as 


s r ({j e }.o) = 


= E ] [ ip- -n 

{m e } e v 


(23) 


with the orientation sign e v e . = —1 resp. 1 recording if the edge is oriented inwards (v = t(e)) resp. outwards 
(v = s(e)). 
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For instance, Wigner’s 6j-symbol can be defined in this way, on the (appropriately oriented) tetrahedral graph, 


J1 J2 J3 
\k 35 je 


= V” 3i-rni / 31 32 3 3 

^ \mi ro 2 m 3 


m 1 , 777 , 2 , 7713 , 
7714,7715,7716 


Jl J5 J6 
-mi -m 5 m 6 . 


J3 J4 J5 
-m 3 -m 4 m 5 , 


32 36 J4 

-m 2 -m 6 7774, 


(24) 

From the definition, it can be seen that flipping the orientation of an edge with color 2 j e changes the spin network 
evaluation by a sign (—l) 2je . Our definition thus really seems to depend on the choice of an orientation o. However, 
we will show in the next section that for a fixed graph and a fixed coloring, the evaluation is constant on the class 
of Kasteleyn orientations (and can be explicitly related to other definitions of spin network evaluations). Then for 
the sake of simplicity, we will drop the o from the notation and simply write s r ({j e }) when there can not be any 
confusion. 

In the context of quantum gravity, the spin network evaluation is interpreted as the scalar product of the spin 
network wave-function with the flat connection state. This actually gives the physical solution for 3d quantum gravity 
(on S 2 x [0,1]) and more generally for topological BF theory. Considering an arbitrary gauge-invariant function ip, 
we can decompose it on the spin network basis states ip^j j upon integrating it against the SU(2) characters: 

V’ttS'e}) = ^j^e) , V’OO = Yl^{3e} sT {{3e}), J Xj* (0e) ^({tfe}) = jj + ^ 1p{j e } ^(O'e}) ■ 


{je} 


{ 3e} 


Due to the orthogonality of the 3j-symbols, the scalar product between two states ip and <p simply reads in this basis 
as: 

1 




{je} 

We now introduce the gauge-invariant flat state, physical solution for the topological BF theory: 


^({ffe}) = / [dh v ] W^ e) g e h t(e) ). 


(25) 


Carefully tracking the inverse group elements and the signs by the following identity on Wigner matrices: 

D^h- 1 ) = (D\h.))l n = = (-1 )<- n - m >DL n ,- m (h ), 

we get the simple decomposition of the flat state and scalar product property: 

fi ({3e}) = ^1°) = (n| (p) = ^s r ({j e }) <p{ je} = (p( I). 

fie} e He} e { ie } 

Spinfoam transition amplitudes between spin network states are then typically constructed from such spin network 
evaluations or projections on the flat state. 


B. Choice of Normalizations and Orientations 


In order to relate our definition of “unitary up to a sign” evaluations to other standard definitions, a couple of 
recoupling identities are required. 


The following equality is equivalent to the definition (24) of 6j-symbols (the equivalence can be proved using the 
properties of 3j-symbols listed in the previous section - see also formula 6.2.6 in [I]), 


E 

77712 


12 — 77112 


Jl J12 J 2 .712 J J 3 

mi mi 2 m 2 ) mi 2 m m 3 


= (-i) 2jl E( 2 ^ +!) r ■ 


3 1 3 2 3 12 


3 3 3 3 23 . 


(- 1 ) 


ji+jz+ja+j 


m 2 3 


3 1 3 
mi m 


3 23 

-m 23 , 


3 2 3 23 3 3 

,m 2 m 23 m 3 . 


(-l) j 


(26) 







11 


Note that the sums on both sides are trivial (the only non-vanishing contributions come from TO 12 = —mi — m 2 and 
TO 23 = mi + m). One can also flip the signs of mi, m 2 , m 3 and/or m and add factors like (—l) J1_mi on both sides of 
the equation without any more changes. We will also need the well known orthogonality relation, 


m 1 ,m.2 


(-1) 


jl — m 1 +j 2 -m 2 +j-m 


(- 1 ) 


j 1+32+j 


J 2 

3i J \ 

Ji 

3 2 

J 

\-m 2 

—mi —m I 

\m 1 

m 2 

777/ I 


3j,j' 

2 j +1 


(27) 


The above equalities can be used to modify the summand in the evaluations (231. Equation (26 1 applies to any 
edge (with spin j 12 in the notation of (26)), while (27) gets rid of faces of degree two. They translate graphically to 



(28) 


(29) 


These equations are unchanged if one flips the orientation of a half-edge on both sides. 


Notice that we have singled out the sign (—l) 2jl on purpose in the right hand side of (26) and (|28|), as this is a 


manifestation of the reason why we call our definition (23) of spin network evaluations “unitary up to a sign”. 

The standard definition of spin network evaluations, which we call the integral normalization and write (T, {j e }) Int , 
is due to Penrose El and is based on associating to each colored graph a linear combination of union of curves lying 
in 5(T) and then evaluating each union of k curves to (—2) fc . Some renormalizations of Penrose’s definition known as 
skein theoretical evaluation and unitary evaluation also exist. As we will prove in Theorem |III.4[ the relation of the 
above “unitary up to a sign” evaluation with the integral spin network evaluation is the following, if o is a Kasteleyn 
orientation on T: 


(r,{je}) Int 


S(r, {je},0) 1 


' ILX-^ + i)! 

FLO*. - 2 je)! 


<r,{je}) sfe ' 


fig 3 a- 

n e (2/e)!‘ 


(30) 


J v denotes the sum of the three spins incident to v. Furthermore, the unitary evaluation is such that the value of any 
theta graph is 1 ; the normalization we use turns out to coincide with the unitary normalization up to a sign which, 
if o is a Kasteleyn orientation on T is computed explicitly: 

s(r,{j e },o) = (r,{j e }) u (-i)i^J», 


Our main object of interest in the present article is the generating function of integral evaluations of spin networks 
with edge variables {Ye} ee p. Provided o is a Kasteleyn orientation so that (30) holds, it can be written 


z Spm { r,{Y e }) = ]T<r,{j e }) Int H Y ? j * = £ 

{ie} e {. j e } 


UeMv - 2je)! 


3 ({ie}, 0 )n F e 2je • 


(31) 


C. About the Tensorial Definition of Unitary Spin Networks 


Observe that if o and o' are orientations on T and o' differs from o on a single edge e then s(P, {j e },o') = 
s(P,{j e },o)(-l) 2 A. In contrast, the standard spin network evaluation does not depend on any orientation of the 
edges of the graph. So the tensorial definition (23) does not coincide in general with the standard one (or it does only 
up to a sign) if the orientation of T is chosen randomly. 

In |12] (section 2.1) a different tensorial definition was provided based on supersymmetric vector spaces (instead 
of even ones) which allowed for a definition independent on the choice of an orientation of the edges of T. Roughly 
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speaking, one stipulates that the vectors of the representation with spin j e have parity 2 j e . Then to each vertex 
one associates the (even) tensor: ij 1 j 2 j 3 (as defined above) and to each edge colored by j e the (even) tensor oj e := 
Y^m e =-j c (~® fi-me ( w h ere <^ e (e)( e ) = 1 if k = m e and 0 else). Observe that because of the parity 
assumption in truth, in order to define uj e we did not use an orientation: indeed applying the flip (exchanging the two 
tensors in w) gives w' = Sm e (~l) 2 j '(~l) je_me ^-m e ® = J 2 m' (-l) 2 je+ 2 m ' +je_m =5^, 0 = u e (where in the 

second equality we just reindexed the sum and observed that 2 j e — 2 m' e is even). Then one defines the spin network 
evaluation as the total supersymmetric contraction of the tensor <^) p 0 (£) v ij*j*j* where, again, in the contraction 

one should take care of introducing signs (— l) ab whenever permuting tensors of degrees a and b. In particular, observe 
that ij 1 j 2 j 3 = (— l) 2jl ( 2 1 2 + 2 1 3 ) i : j 2 j 3 j- l but this apparent asymmetry is compensated during the contractions. 

In our case we deal with the orientation problem by choosing a Kasteleyn orientation on T (which is supposed to 
be planar). Observe that by the admissibility condition of the spins (i.e. the fact that for each vertex v the sum J v 
of all the spins of the edges surrounding v is integer) and by the above observation on the behavior under the switch 
of the orientation of an edge, formula (231 provides a well defined function on the space of orientation classes which 
are defined as follows: 


Definition III.3 (Orientation class). Two orientations on T are equivalent if they can be obtained from one another 
by a finite sequence of moves consisting in switching all the orientations of the edges incident to a vertex. We will 
denote the equivalence class of an orientation o by [o]. 

In our case T is supposed to be planar and we can equip it with a Kasteleyn orientation; there may be more than 
one such orientation but it is easy to see that all of them are equivalent in the sense of Definition |III. 3| hence the spin 
network evaluation on T is well defined. 


Theorem III.4. Let T C R 2 be a planar connected trivalent graph and o be an orientation on T. Then, the following 
holds for each coloring {j e } of the edges of T if and only if o is a Kasteleyn orientation: 


(T,{je}) Int 


+ 1 )! 

riei/'^ — 2je)! 


S( r ,{je}, 0 ). 


Proof. We will first prove that among all the orientation classes on T there is at most one, the Kasteleyn orientation 
class, which can have the required property. The idea is to compare Equation (30) with the standard evaluations 
only on “curve colorings” i.e. colorings with values in { 0 , ^}; these colorings are naturally identified with the set of 
curves embedded in T (the curve associated to a coloring being the set of edges with spin f). Recall that the standard 


evaluation of such a coloring is (— 2 ) fc , where k is the number of connected components of the curve associated to the 
coloring. 

The proof is actually similar to the (constructive) proof of the existence of Kasteleyn orientations. Pick up a 
spanning tree T C T, whose orientation is the restriction of o to the edges of T. Then we show that requiring the 
curve coloring evaluation to be —2 around each face of T enforces an orientation on the edges of T \ T, and it is such 
that one gets a Kasteleyn orientation on T. 

We denote T* the dual graph to T, whose vertices represent the faces of T and an edge connects two vertices dual 
to the faces /i, whenever there is an edge of T which belongs to the boundary of both faces. Therefore an edge of 
T* uniquely identifies an edge in T and the other way around. We denote (T \ T)* the dual to T \ T. It is a tree: it 
is connected because T has no cycle and it has no cycle since T is connected. Since it goes through every face of T, 
(T \ T)* is a spanning tree of T*. We choose a root vertex (for instance dual to the outerface) and consider a leaf /*, 
i.e. a vertex of degree 1. It is dual to a face / whose boundary edges all belong in T but one, say e/ (which is dual 
to the edge incident to the leaf). We denote the boundary edges of / by 


df = {df H T) U e f , 


which corresponds to the partition of the edges T = (TnT)U(r\T) restricted to df. We want to check that imposing 
s(df,o) = —2 implies for e/ an orientation such that there is an odd number of clockwise edges in df. By switching 
the orientations as in Definition |III.3| we may assume that all edges of df D T are counter-clockwise around /. 

Let us start observing that 


*n° - ^ 2 ( i 


1 en - e* 


3 ), and 


1 1 


- 0 - 

o o 


— 7 =(e? ®eS®e 2 ! 

V 2 2 u -2 



0 e° 0 e 


1 



2 


then if n is the number of vertices (or edges) contained in df then the evaluation is 


1 

7 ^ 


n 


(i + i)- 


s{df,o) = ± 
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FIG. 2. An example of how to reduce T to a connected sum of theta-graphs (n = 4). 


Indeed there are exactly two non-zero summands in Formula (23) and they can easily be checked to have the same 


sign. Since the normalization factor )! equals \/2 n , it is clear that we get —2 exactly for one of the two 

possible orientations on e/ and 2 for the other. Then a direct inspection in Formula (23) and in the tensors iii 0 and 
shows that 


Uni 

o u o 


s(df,o) = 2 sign(e/) (-1)“( 9/ ) 

where sign(e/) is 1 if ef is oriented pos itive ly by df and —1 else and a{df) is the number of “large angles” 
as defined in Lemma 


II.4 


By Lemma II.4 since df contains no internal vertices, we have (— = 1 ; 


in Of 
hence 

s(df,o) = —2 if and only if e/ is oriented negatively with respect to /. Switching back the orientations on df fl T 
to o in the sense of Definition |III.3[ we preserve an odd number of clockwise edges, as claimed. A straightforward 
induction from the leaves of (T \ T)* to its root concludes this part: only the Kasteleyn class can provide the standard 
evaluation for all colorings of T. 

To prove that it actually does, for every coloring {j e } on T, we start by observing that the claim is true if F is 
the theta-graph. Moreover if T is generic and equipped with a Kasteleyn orientation o and a coloring {j e }, then we 
can transform it into a graph which is formed by a connected sum of theta-graphs. This is done by means of a finite 
sequence of Whitehead moves similar to Equation (261. To see this one argues by induction on the number of vertices 
of F: choose a region of R 2 \ F with, say, n edges and pick one of them up, say e. Applying n — 2 Whitehead moves 
one can “slide one endpoint of e” along the boundary of the face until one gets a new graph containing two arcs (one 
of which is the image of e under the sequence of moves) having the same endpoints. Then the so-obtained graph is a 
connected sum of a theta-graph with a “simpler” one. Figure [2] is an exemplification of the case n = 4 of this process. 


The algebraic translation of a Whitehead move in terms of spin networks, for our unitary up to a sign definition, 
is precisely Equation (28). Observe that indeed, one can always orient the edge carrying j ’12 as in the left hand side 
(l.h.s.) of (28), by switching the orientations of the edges incident to one of its endpoints, if necessary. However, the 
Whitehead moves does not preserve the Kasteleyn orientation: the drawing on the right hand side (r.h.s.) of (28) is 
not equipped with a Kasteleyn orientation anymore. The face going along (j,ji 2 ,ji) on the l.h.s. loses a clockwise 
edge in the process, while the face along gains one. Therefore this can be corrected by switching the orientation 

of the edge carrying j\ on the r.h.s. This brings up an extra factor (—l) 2jl which cancels the one we already had in 
(28). The Whitehead move respecting the Kasteleyn orientation class is thus 



In this equation, the orientation of the half-edges can be flipped, as long as it is done on both sides. 

Similarly observe that if o in the left hand side of Equation (29) is a Kasteleyn orientation then o', on the r.h.s. is 
not. To correct this, we need to switch the orientation of the only edge (colored by j) thus introducing an additional 
factor (—1) 2 F 

Now the crucial observation is that the standard integral evaluation satisfies the same equations, up to the fac¬ 
tors A(a, b, c ) = (a+b _ c ;^^|- + c-a)! which we have not taken into account yet. Indeed letting A(a, 6, c) re¬ 


normalize s(r, {j e },o) to the standard normalization of (T, {j e }) Int we have to multiply the left hand side by 
\J A (ji ■ j 2 , ji 2 ) A(j '3 , j, 2 ) and the r.h.s. by yj. A (j 1 , j, j 2 3 )A(j 2 , h ,J 23 )• But, taking into account Racah’s for¬ 
mula for Wigner 6j-symbols and the fact that the standard evaluation of a theta-graph colored say by j\,j,jis is 
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(—l)J' 1 +J+i 23 A( i 7 i, j, J 23 ) we see that the standard evaluation (T,{j e }) Int satisfies the following recoupling identities: 



I 31 

V A U2,j3,j23)A(ji,j23,j) 



(33) 



(-l)n+n+i(-l) 2 i6 j j,A(j,ji,j2) 

2 j + l 


(34) 


Since the values of vn V (J V + l)!/ru<4 — 2j e )!s(r, {j’e}, o) with o a Kasteleyn orientation satisfy the same re¬ 
coupling identities as the values of the standard evaluation (r, {j e }) /nt , the theorem is proved. III. 4 


D. Gaussian Representation of the Generating Function and its Loop Expansion 


Similarly to the Ising model, we will provide a representation of this generating function for spin network evaluations 
as a Gaussian integral, but this time in terms of usual even-Grassmannian integration variables. 

Let r be a planar graph, o a Kasteleyn orientation on T and 


Z s P m (T,{Y e }) = Y / 

Ge} 


' IL(’A. + 1)! 

~ 2 j e )! 


IIn 2j ” = E( r -ije}) Jnt U Y e je 

e be} e 


(35) 


By Theorem m ~A\ we can omit the notation o. In this section we provide a Gaussian integral whose value is 

Z S Pin(T,{Y e }). 

A similar integral was provided first in m The difference between the approach in the present paper and that 
in [121 is that here we use complex-valued integrals and we deal with a Gaussian integral which is convergent as it 
is, while in [2] a suitable regularization procedure was used to make sense of the integrals. On the downside the 
integral considered here is on a space whose dimension is twice that of the space used in [T2]. Another Gaussian 
integral expressions was given in sum, but it involves integrating variables living on the vertices of T instead of 
its half-edges as done here. The relation between these two formulas and their equivalence will be examined later in 
Section IVII 

One of the interesting points of having such Gaussian integration formulas is that they allow to re-prove Westbury’s 
theorem m asserting that for a planar graph T, 


Z s P in (r,{Y e }) = Pf 2 where P r = £ ]jY e , 

c curves in T eEc 


where “curve embedded in T” means a (possibly empty) union of edges of T homeomorphic to a union of circles. It is 
worth noting that Westbury’s initial proof was based on a totally different approach, namely computing spin networks 
via chromatic evaluations. 

Let us first fix some notation. Recall that for an angle a £ A of a planar graph T, we let s(a) and t(a) be the half 
edges at the source and target of a (using the cyclic ordering of the edges around the vertices) and X a = -^/Y s ^Y t ^ a y 
For each half-edge h contained in an edge e and incident to a vertex v, let z ev ,w ev be complex variables and z ev ,w ev 
be their complex conjugates. 

Theorem III. 5. With the above notation it holds: 

Z Spm {T,{Y e })= [ 

J(C 2 ) 2 # E 


n 


dz P ,,dz P ,,dw P ,,dw, 


ev e -£ E „(l 2 


v\ 2 + \Wev\ 2 ) p - E I (2.(e)»t( t )-'».(,)2i(e))+E a (z o(ct) ( o) -1U o(a) Z t(a )) 


where the r.h.s. is a well defined integral on <CA# E . Computing the value of the integral one 

p - 2 

- t r • 


also gets Z s P in (T, {Y e }) = 


















15 


Proof. Let us prove the first statement. We start with the generating function of Wigner 3j-symbols, which is known 
[34] to be given by: 


3 1 02 03 
m i ,7722 

where ./ = jj + j 2 + js and a denotes an angle. In the notation if a is the angle between the edges i and i + 1, then 
Z S {a)W t {a) - Ws(a)Z t (a) = ZiW i+1 - W^ i+ 1 , for * = 1,2, 3. 

Taking the product over all the vertices of T of the above generating sums we do not yet get the desired generating 
sum for two reasons: 



Jl J2 J3 
J7ii ro 2 m 3 






i=1 \/( j - 2 ji)Kji - miV-Ui + m,)! 


= exp^X a (0 5(ct )M; t(a ) - w s ( a )Z t ( a )). (36) 


1. We need to consider only products where the spins and magnetic numbers of half-edges contained in the same 
edge match; 

2. Even after solving point [T] the total weight of the state associated to a spin network with a choice of magnetic 
numbers per each edge would be off by a factor J][ e (—l)- 7e_me {{j e — rn e )\{j e + w e )!) _1 (compare with Formula 

We solve both preceding problems by considering for each edge e € E the following (well-defined) Gaussian integral: 


dztdztdz s dz s dwtdwtdw s dw s , 


z t \ 2 + \z B \ 2 +\w t \ 2 + \w B \ 2 ) -(w t z B -z t w B ) Jt+mt j B -m a j B +m a jt~m t 


J C 4 7r 

= dj u j B d mt -m B (jt + m t )\(jt - m t )\ 


(37) 


Thus taking the integral over all C 4 # E of the product over all the vertices of the generating series (36) and the 
product over all edges of the above Gaussian factors we get the desired result. Observe that the so-obtained Gaussian 
integral is convergent when {Y e } are near 0. 


One could then try to switch back from complex variables to real variables, and avoid the doubling of the number of 
integrals, by using the Segal-Bargmann transform, which maps complex holomorphic monomial z n onto the Hermite 
polynomials H n {x). This works and leads to a real Gaussian integral but finally produces a non-trivial and non-linear 
action in terms of the couplings X a . 

We now prove the last statement. The quadratic form associated to our complex Gaussian, is given by a (8#E x 8#E) 
square matrix. In terms of the real parts of the z and w variables and their imaginary parts it reads as: 



—A T +B r + 21 

M(A r + B r ) 

~i{A T + B t ) 

A T ~B V + 21 


where A r and B r are the symmetric real matrices of size 4#E defined as follows: 

ZA Z — 2 ^ ' l ^s(e)Zt{e)^)i ZB Z = 2 ^ ^ X a (z s ( a ')Wt( a 'j ' l ^s(a)Zt(a)') ■ (^) 

e a 

Using the facts that (4 r ) 2 = I and 4 r U r = — B r A r (which are easily checked by the definitions of A r and B r ), 
and summing first — * times the last columns to the first ones and then i-times the new last lines to the first ones we 
compute: 


detM r 
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2i(I - B r ) 

—2i(I + A r ) 

A r ~B r + 21 


4 4#e det(I + A r B r ) = 2 8#e det(4 r + B r ) 


where the second equality holds from the fact that the determinant of a 2 x 2 -block matrix (^ ^) is det(M 4 — 
M 2 ) det(Mi) if M\ is invertible. 

To conclude we now follow the lines of the approach used in [12] : observe that in terms of the real parts of the z 
variables and of the real parts of the w variables both A r and B r are symmetric and of the form A = ( ^cT ) anc ^ 

B = ( 0 , B 0 * with A 11 = —A and B n = —B. Thus det(4 r + B T ) = Pf [A’ + B') 4 where Pf denotes the Pfafhan. We 

claim that Pf {A' + B') = Pp- To interpret it we use Kasteleyn’s approach 35] by interpreting the matrix A' + B' as 
the adjacency matrix of the oriented graph T 7 obtained by replacing each vertex of T with a triangle and orienting the 












16 


new edges counter-clockwise, while keeping the edges coming from those of T oriented as in R It is straightforward 
to check that R inherits then a Kasteleyn orientation. 

Remark that the edges of R are of two types: those corresponding to edges of R say of “edge type” and weighted 1, 
and those corresponding to angles of T, say of “angle type” and weighted X a . By Kasteleyn’s results the value of the 
Pfaffian is the generating sum of the dimer configurations on R weighted with the weights of the edges contained in 
the dimers. Dimer configurations on R correspond bijectively to curves in Ik consider the angle-type edges contained 
in the dimer configurations; they can be completed in a unique way to a curve in R by adding some edges of edge 
type; then retracting the triangles of R to points we get the desired curve. Reciprocally each curve in T meets a 
vertex either 0 times or 2 times: in the latter case it identifies an edge of angle type in R, in the former we associate 
to it the three edges of edge type of R surrounding the triangle corresponding to the vertex. 

Since the orientation of R is a Kasteleyn orientation, we conclude by the results of [35] that Pf {A! + B') is the 
generating series of dimer configurations on R weighted as above, which in turn equals ^ C crriacc^“- Recalling 


This theorem implies the following remark, central to our present work, on the duality between the 2D Ising model 
and spin network evaluations: 

Remark III. 6. Let T be a planar, connected and 3-valent graph. Then if we match the Ising couplings y e to the 
parameters of the spin network generating function by Y e = tanh y e , the following equalities hold exactly: 


that X a = \/Y s / a )Yt(a)> we are done. 


III.5 


Zf(T,{X a }) 2 Z Spin (T,Y e ) = 1, z Ising (T, {y e }) 2 Z Spin (T, {Y e }) = (^2 2#V Q cosh(y e ) 2 J , (39) 

where we keep the matching X a = v/R s (a)^t(a) between angle and edge variables ensuring that Yl a &c ^ a = n eec Y * 
around any loop C in the graph T. 


IV. DUALITY THROUGH SUPERSYMMETRY 


In order to understand the relation between the Ising model and spin network evaluations, one can play a little game. 
On the one hand, the Ising partition function on a graph T has a simple loop expansion as Pr = 1 + Py = 1 + X] 7 cr W 7 
(up to pre-factors) where we sum over all non-empty even subgraphs 7 . These even subgraphs are simply identified 
as unions of disjoint loops when T is 3-valent. The amplitudes W 7 are the weights J ([ eg7 ^e- On the other hand, the 
generating function of the spin network evaluation is related to the squared inverse of P 7 , which can be expanded as 
a power series: 


P 


-2 


1 

(1 + Pr ) 2 


£(-!)"»/£ = £(-!)"« E W 71 ..W 7 „. 

n£N n£N 7i--7n 


(40) 


Such union of even subgraphs can be considered as an arbitrary subgraph G with extra integer labels c e on each edge 
indicating how many times an edge e belongs to one of the even subgraphs 7 ,. We have the obvious constraint that 
the sum of c e ’s around each vertex v is even, plus the less obvious constraints that they must satisfy the triangular 
inequalities. These colors c e are actually to be thought of as twice the spins j e . This allows to write the power series 
above as: 


P r 2 = l+E E C G (c e )W^ } , 

GCr {CefeEG 


Y'r 1 


nv 

eGG 


(41) 


where the coefficients C G (c e ) count, up to the pre-factor (—1 )™n, the number of ways we can decompose the subgraph 
G with admissible coloring {c e } as a union of even subgraphs 7 j. This is actually the expansion in terms of spin 
network evaluations. But it explains, in some sense, how while the Ising model sums over configurations with the 
occupation number on its edge being 0 or 1 , and can be thought as fermionic, the inverse of its partition function 
sums over configurations with arbitrary integral occupation numbers (with some non-trivial combinatorial factors, 
which are actually the spin network evaluations) and can be thought as bosonic. We formulate in a rigorous way 
below defining a supersymmetry relating the Ising partition function and the generating function for spin network 
evaluations. 
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A. A Supersymmetric Theory 


We showed in Proposition II.3 that the Ising partition function (on the graph P) can be formulated as a Gaussian 
fermionic integral. It directly evaluates the determinant of the quadratic form, which is simply the polynomial Pp. 
On the other side, the generating function of spin network evaluations is, by Westbury theorem, the inverse squared 
1/P 2 . It is thus given by a bosonic Gaussian integral to get the inverse determinant, together with a doubling of the 
variables and an anti-symmetrization, in order to reach the correct power of the determinant. This is exactly what 
we have achieved in the previous section. 

Let us now start with the following “meta-theory”, with both fermionic and bosonic degrees of freedom, with a 
copy of M 2 ® (/\ R ) 2 on each half-edge, Z(T, {Y e }) = Zf(T, {y e }) 2 Z C (F, Y e ): 


z(r,{F e }) = / 5 n 


half-edges h 


dz h dw h d4’ h dlj h ex P ^(z s ( e )W t ( e ) - W s ( e) Z t ( e) ) + + Vs(e)Vt(e 


exp 


x>-( 




v)Zt(a) “1“ '^Ps(a)'4^t(a) “1“ Vs(ct)Vt(<y) j I (42) 


with X a = JY s[a) Y t ( a) as before. The variables z, w are bosonic and are the spin network degrees of freedom, while 
the odd-Grassmaniann variables ? 7 , 1 / are fermionic and are the Ising degrees of freedom. Obviously the two families 
of integral factorize. One recognizes the Grassmannian integral of Proposition II.3 doubled (integration over 1 f> and 
rf). As for the bosonic part, it is a real version of the formula of Theorem III.5 which is however not well-defined 


(its quadratic form has negative eigenvalues). This Gaussian integral should in principle be constant, thus leading 
to the result (39). Although the above integral is ill-defined, it is a good starting point to define the supersymmetry 


transformations and later extend them to the well-defined complexified action. 

Notice that the argument of the exponential splits into contributions labeled by pairs of adjacent half-edges i,j , 
either on the same edge e or on the same angle a; in the former case we suppose i = s(e),j = t(e) and in the latter 
i = s(a),j = t{a). The total action thus reads S = Y/ 
are adjacent, and 


<i j > A ijSij where < i,j > means that the half-edges i and j 


Sij - (ZiWj - UHZj) + ipii>j + ViVj , (43) 

with the couplings X^ = 1 on an edge and Xij = X a as before for an angle. 

We define the supersymmetric operator Q of odd parity (i.e. an odd derivation) acting on every half-edge i on 
Si = R 2 ® (AR)? =R - 12 as 

Qzi = ipi, Qnh = Tji , Qipi = Wi, Qrji = -z t . (44) 


We naturally extend the action of Q to the algebra of functions on ®jSi via the graded Leibniz rule. It then follows 
that for any adjacent half-edge j , we have: 


QSij = ipiWj + Zirjj - 4>jWi - Zjrji + Witpj - 4’i w j - ZiVj + Vi z j = 0, (45) 

thus the action S = Yl<h g > XhgShg is Q-closed for each half-edge i. 

Moreover, the action is actually Q-exact (but Q 2 does not vanish), 


Sij = Q$ij = with = Zi4)j + WiTjj. (46) 

Viewing the integral Z(T) = f exp (J2 <h g> X hg Sh g ) as a function of the couplings {X^} associated to pairs of 
adjacent half-edges, we write 


dI(T) 

dXij 


/ Wdzhdwhd'tphdr/h exp( Y] X hg S hg ] = / T\dz h dw h di/j h dr]h Ql $ij exp( X hg S hg )). (47) 
•. K <h, g > J -Is h V V, 


<h,g> 


Using the lemma |IV.1| proven below, stating that the integral of any Q-exact function vanishes, we can deduce that 


dZ(T) 

dXij 


= 0, 


(48) 


This proves without explicit evaluation that the integral Z(T) is a constant. Since it is independent of the variables 
{X^}, one can then evaluate it on any set of values making the integral simple. 
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Lemma IV.1. Let f £ C°°(K^ 2 ) be a function admitting a Berezin integral. Then it holds: 

/ dzhdw h difhdr\ h Qf = 0 . 


(49) 


Proof. Let e be an odd-Grassmannian variable, e 2 = 0 and perform the parity-preserving change of variables: z K > 
(1 + eQ)z, and similarly on Ordering the variables as z,w,ip,r], the Jacobian matrix of the transformation 

reads as: 


J = 


(1 o 

0 —e \ 

0 1 

e 0 

e 0 

1 0 

\0 e 

0 1 ) 


Its Berezinian (or super determinant) is simply evaluated 


SdetJ = det 


I - 


0 -e 
e 0 


e 0 
0 e 


= 1. 


Applying this change of variables on the integral of / gives: 

/ TT dz h dw h dif h dr) h f = / TT dz h dw h dif h dr] h ( / + ^ c(Qz h d Zh f + Qw h d Wh f + Q4’ h d^ h f + Qvhd Vh f ) ) ■ (50) 
J S h JS h V h ) 

To get this equation, one has to be careful with the parities of e,Qz,Qw and similarly that the operators djj,,d v are 
graded. Then observing that 


Qz h d Zh f + Qw h d u ,J + Qifhd^bf + Qvhd Vh f = Q(f), 


the desired result follows. 


(51) 

TvT 


We now would like to follow this procedure for the well-defined complex Gaussian integral for the generating function 
of the spin network evaluation given in Theorem EH To allow to match bosonic and fermionic variables, we also 
use the formulation of the Ising model in terms of complex fermions and we then use Formula (15). 

Thus we consider the partition function: 


z c (r,{i;}) = 


r 4#E 


/ TT dz h dzhdwhdw h dtphdi> h drihdr] h e +KI ) 

^ half-edges h 


E e ,„ (^eVl+^eVe) “Ee (vb(e)&(e) +Vs(.e) Vt (e) ) W (z„( Q ) W t (<*) - Zt(c) +bs( Q ) bt(c) +Us(, a ) ) 


, (52) 


We consider the same supersymmetry generator Q as before and we extend its action on the complex variables, both 
bosonic and fermionic by assuming its compatibility with the complex conjugation, Qf = Qf. Then we decompose 
the complex path integral distinguishing the half-edge amplitude, the edge action and the angle action: 


Z c (r,{A e ,„, M e,^}) = 


r 4#E 


][| dz h dz h dwhdw h d'ifhdi> h d’i]hdr] h e 


'n. P Ee Se S a 


(53) 


(54) 


° half-edges h 

where we have added the extra couplings A ei „ and [i e , which can be set to 1 , and 

K e , v = KI 2 + |<| 2 -«-«. 

= ^s(e)^t(e) ^s(e)^t(e) H” '*Ps(e)'$t(e) T Vs(e)Vt(e)i 

Sql = Zs(a.)^t(a) ^s(ot) (a) T T ^?s(o:)^ 7 f(Q') • 

A quick calculation shows that all three types of terms are invariant under supersymmetry. They are Q-closed and Q 
exact: 


QK etV = QS e = QS a = 0, 


-^e,v — Q i^Pe,v^e,v 

= Q (^ z s(e)'^t(e) T ^s(e)^t(e)^ 
$ot Q T VJs(a.)Vt(a )) 


(55) 
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Using Lemma IV.1 the integral (53) does not depend on any of the coupling constants A, p or X: 


dZ c (T) _ dZ c (T) _ dZ c (T) 

&Pe 


For instance, we can set A = p = 1 as in the original integral and then set all X a = 0 (or equivalently all Y e = 0). 
This shows that the path integral is equal to the case where the angle interactions are killed and we are left with a 
free system with decoupled edges. We can also send the couplings A ej „ to +oo. This is the idea of localization: the 
integral localizes on z = w = 0 and ip = p = 0 and has a trivial evaluation. 

This opens an interesting direction for future investigation, on how this localization due to the supersymmetry 
constrains the value of the observables and correlation functions of both the Ising model and the spin network 
evaluation. 


B. Towards Non-trivial Coupled Supersymmetric Theories 


Of course, the above integrals representing the Ising partition function and the spin network generating function are 
Gaussian, and therefore localization is not really necessary since we know how to compute these integrals explicitly. 
We can nevertheless push the logic further and propose non-linear extensions, beyond the quadratic action, which 
would still be supersymmetric and thus localizable. This provides a rather large class of non-Gaussian integrals, which 
we can now compute exactly using our new tools, coupling non-trivially the spin network evaluations and the 2D Ising 
models on a planar graph. 

Since each terms of the action, K e:V , S e and S a are Q-closed and Q-exact, so are arbitrary powers of these terms. 
This means that we can add arbitrary powers of each of them to the action, keeping the integral invariant under the 
same supersymmetry as above and therefore localizable. This leads to a whole class of explicitly computable integrals: 


Z c 

n,p,q 


(r, {A ei „ 


ftej W}) 



A e ,„if e , 


-E e ^eS e +E a -Y c> S Qe -E e ,„ 


X i P ,l K e,v~T, c + 


(56) 


where n,p,q are arbitrary integers, larger than or equal to 2, and the measure dp is the product measure over all 
half-edges h of all the even and odd Grassmann variables dp = dz] l d'Zhdwhdwhd'iphd'ip h dphdrji l - 

As before these new integrals are constant and do not depend on the specific couplings pi^ and X^' 1 . In the 
original Gaussian integral without the new higher order terms, the half-edge terms K e v were the Gaussian measure 
factors, the edge terms S e were the anti-holomorphic part of the action, while the angle terms S a were the holomorphic 
part of the action. In a sense, we were taking the scalar product between the edge action and the angle action, thus 
gluing the angles with the edges using the Gaussian measure. In this context, it does not seem useful to add non¬ 
linear terms to all three part of the action and we focus on adding higher order terms to, say, the angle action, thus 
considering: 


2 £(r,{A ei *,Me,n}) 



\ e ,vKe,v—^2 e / J 'e5'e+2^ Q , X a S a 




(57) 


where we can compute the arbitrary power S™ for n > 2 taking into account that the ip' 1 s and ?y’s are odd- 
Grassmannians: 


S 


n 

a 


^s(a)^t(a) '4^s(cn)'4^t(a) 

^s(ot) ^t(a )) ”1“ ^ ^s(a)-^<:(a)) {.'4^s{a)'4^t(a) “I - Vs(ot)'nt(a.)') 

nilTj l) / \ n —2 /it \ 

S (^s(a)^i(a) ^s(a)^t(cx)j lr's(a)ri(Q)^s(a)^(a)) ■ 


(58) 


The first term corrects the weights of the generating function of the spin network evaluation. This will modify its 
stationary points, as explored in Section [VI B[ and thus the geometrical background. The second term then produces 
a geometry-dependent coupling for the Ising models, while the third term creates a coupling between the two Ising 
models. 

Let us focus on the quartic case, n = 2. Then the coupling between Ising models is not geometry-dependent. We 
can write this integral explicitly: 


-Z?(r,{A e>v ,p e ,Y e }) = J s dpe 


~ ^ ~l e v ^e,vKe,v — ^2 e Me'S'e+^o: X a S a 


(59) 


Ta4 2, [('.(»)®l(a)-».W2t(«)) 2 +2(2,( n| »),( 0 )-1«,(a | ^(a))('f«(,)^( n )+'l,W%w)+f,(«)*(«)'l.(«)%(»)] 
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( 2 ) 

Despite the new coupling X„ , we know that this integral still evaluates to 1 as before. It would be enlightening to 
investigate how the generating function for spin networks is modified and what would be the physical effect of the new 
couplings for the Ising models, which seem to allow us to study the Ising models with dynamical, geometry-dependent 
couplings. 


V. THE INTERPLAY BETWEEN SPIN NETWORKS AND ISING CORRELATIONS 

In this section, we investigate the link induced by the duality developed above between the correlations of the 2D 
Ising model and the probability distribution defined by the spin network evaluations. 


A. Coherent Spin Network States, Geometric Interpretation and Criticality 


In the context of quantum gravity and quantum geometry, following the logic developed in mi. coherent spin 
network states with good semi-classical properties can be interpreted and effectively seen as generating functions 
for spin networks in the spin basis. It is similar to the way coherent states for the harmonic oscillator define the 
exponential generating function for the eigenvectors of the number of quanta operator. 

We introduce the wave-functions on the graph T depending on group elements g e £ SU(2) living on the edges e £ E: 



/ n- 


(60) 


where the use the following convention for spinors | z) and their duals \z], as introduced in [36l 137] : 


l*> = 



£ C 2 , 




and d A z := dzdzdwdw. 


This allows to encode both complex variables z and w in a single complex 2-vector |;s). The SU(2) group acts naturally 
on C 2 via its fundamental representation, and the dual spinor transforms exactly the same way as the original one: 


g : \z) ^ g > \z) = g \z) 



1.9 "] =9 \z\- 


The scalar products between states and with their dual provide two bilinear forms invariant under the action of SU(2): 


(2l | Z 2 ) = (Z!Z 2 - WiW 2 ) , [2l|22) = (21 W 2 - 22 ^ 1 ) . 

This implies that the wave-functions ^{g e ) introduced above are gauge-invariant as expected under SU(2) trans¬ 

formations acting at the vertices of the graph: 

V{M„ 6 v e SU(2)# V , <t> c { f a} (g e ) = <t>{x a }(Ki)9ehtu) • 

This set of coherent states is specially interesting for our perspective because their evaluation at the identity 
g e = I, Ve, or equivalently its projection on the flat connection state (251 gives the generating function of spin 
network evaluations: 

= W?l}> = Z s * n {{Y e }) = / nj^ + l )! ({ Je} )]-[y e ^ , 

Del » ZJ e ). 

still assuming the matching X a = \jY s (a)Yt( a )- One can further get the whole decomposition of the coherent states 
<j in the spin network basis j defined in (22): 


= E 

lie} 


I ri L .(^ + i)! 

Wev^ v ~ 2j e )! 




} ■ 


( 61 ) 
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For here, we can consider two types of probability distribution and averages for the spins {j e }- We can look at the 
averages defined by the amplitude 4 >Tx y \ (I)■ that is the averages weighted by the spin network evaluations: 




’) = 


Z s P in ({Y e }) 


{>} 


Y\ev(Jv ~ 2j e )! 


<{je}) H Y e 2j * 


(62) 


This is not strictly speaking a true mean value and a probability distribution since the spin network evaluation s({j e }) 
can be negative. Rather it should be interpreted as an operator insertion in the projection over the flat state: 


or-j'Sr) = 


TO -j 1 


\r oh 


{x~Y 


(: n\(j) coh 


{Xct}' 


(63) 


where the operators j acts by multiplication by the spin j in the spin network basis. These averages will be directly 
related to the Ising correlations, as we will investigate in more details in the next section. They are functions of the 
couplings, Y e which can be thought of as a background geometry in which quantum fluctuations of the spins j ei take 
place. 

On the other hand, we can consider the true expectation values for the spin operators in the quantum state ,): 


(jel-Je^coh = 


<^1} I 


1*11 1 
Jei “Je 


1^1}) 


E e} «- 


■Je^ 


n 


ii c (j»n)i n E Je 

(Jv 2 j e )! lie (2j e +l) 


^{ 1 } 


t^l}) 


E 




fie} n e „(Jv-2je)! He (2j e + l) 


Y? 


(64) 


In practice, this is the expectation of the product of half-integers j ” 1 • • • with respect to the following probability 
distribution: 


P({je}) 


E {jB} P({je}V 


with 


e({«> = n n 

1 1 


( J v + 1 )! 


n 


Y* je 


j{Jv 2 Je)! (2j e + 1) 


(65) 


Notice that it does not feature the spin network evaluation. Let us investigate in more details the shape of this 
probability profile. It will shed light on the properties of the chosen weights for the generating function of spin 
network evaluations and on the geometrical interpretation of these coherent states. To this purpose, we will compute 
a limit by letting jW : E -> N be a sequence of colorings whose spins go to infinity linearly in n, i.e. such that there 

exists a real valued coloring j°° : E —> N such that limn-^ E (oo E =1, Ve £ E. For simplicity, in the following 

computations we will suppress the indices ^ and and study the dominant regime for these weights for large 
values of n using the Stirling approximation for the factorials: 


p({je}) ~ ri 


n 


\j j v +1 


(2j e + 1) 27 t YleBv '/Jv 2y 


n^n 


(Jv + 1) (J ” + 1) 

UeBv(Jv-2 


where we have separated the exponential contribution from the polynomial pre-factors and by ss we mean that the 
limit of the ratio of the two sides when n tends to oo is 1. Assuming that the exponential factor controls the main 
behavior of this distribution, similarly to a Poisson distribution, we focus on the exponent: 


$ = E 4 ^e 111 E + E 


(Jv + 1 ) ln(J„ + 1 ) - E (Jv ~ 2 Je) ln ( J v ~ 2 ie) 


e3v 


( 66 ) 


Let us do a stationary phase approximation, leading to a Gaussian approximation for the distribution p peaked 
around maxima given by the stationary points of <f>. Still assuming that the j’s are large and neglecting sub-leading 
contributions, we get: 

_ Q ^ yA ^ ( Js{e ) ~ 2 j ei )(Js(e) ~ 2 Je 2 ) ( Jt(e ) — 2 jei)(Jt(e ) — 2 Jg 2 ) 

&je Js(e)(Js{e) 2 J e ) Jt{e)(Jt{e ) 2 J e ) 

where ei,e 2 are the edges touching e in s(e) and ei,e 2 are those touching it in f(e), as depicted in Figure [3j At 
this point let us do a little triangle geometry. Assuming that we have three edges of length l , li and I 2 , with total 
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FIG. 3. The notation used on dual triangles to T. 



half-perimeter L , we have Heron formula for the area, a similar formula for the inner circle radius and expressions for 
the sine and cosine of the opposite angle 7 to l: 


A 2 = L(L - l)(L — h)(L — h), r 2 


(L-l)(L-h)(L-h) 

L 


1 + cos 7 


2 L{L-l) 

hh 


sin 7 = 


2 A 
hfa 


Combining all this, in order to describe the stationary point, it seems natural to define dual triangles to the (planar) 
graph F by considering triangles around every vertex with edge lengths given by l e = 2j e (cf. Figure [3]). Then the 
right hand side of © can be re-interpreted in terms of the opposite angles to the edge dual to e in both triangles 
dual to the source and target vertices s(e) and t(e). This gives the critical values of the couplings: 


Y 


( sin 7 | (e) \ 

\1 + COS 7 J 


f sin 7 e^ e) \ 

\ 1 + COS 7e (e) J 



tan 



( 68 ) 


They correspond to a Euclidean structure on the plane obtained by gluing Euclidean triangles dual to the edges of T 
whose edge lengths are 2 j e , up to an arbitrary global rescaling. 

We will not perform here the details of the resulting stationary phase approximation, but we will focus on the 
meaning of these “geometric” couplings from the point of view of the Ising model. In the special case of a honeycomb 
lattice (regular hexagonal), all the angles are 7 = f and this value of the geometric coupling matches exactly the 
value of the critical value of the Ising model, Y c = ^=. 

This actually holds for the much larger class of isoradial graphs, i.e. such that all faces are inscribable in a circle of 
a given radius and whose center is in the face. Let us consider such an isoradial embedding of our planar graph T. The 
dual vertex of a face, i.e the corresponding point of the triangulation, is chosen as the center of its circumscribed circle, 
as illustrated on Figure [4] Each edge of the graph F is at the intersection of two circles, while the dual triangulation 
edge links the centers of those two circles and is orthogonal to the graph edge. In this special setting, it turns out 
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that the half-rhombus angle 9 e associated to the graph edge e (or equivalently its dual triangulation edge) is equal to 
both opposite triangle angles y e = 7e*' e ' ) = 7e^ (which are then equal to each other). Then the critical Ising couplings 
read (see [38],[29] and the survey [30] for additional material): 


e 2 vl = 


1 + sin 9 e 
cos 9 e 


Y e c = tanh y c e = 


e 2y e —1 1 + sin 9 e — cos ( 


e 2y i + 1 1 + sin 9 e + cos l 


= tan y = Y e geom . 


(69) 


This provides a neat geometrical interpretation of both the generating function of spin network evaluations and of 
the critical regime of the 2D Ising model, at least in the context of isoradial graphs and their dual triangulation. It 
would be very interesting to investigate if this correspondence between stationary points in the spin networks and 
critical couplings of the Ising model is more general than this setting and can be generalized by our formula (68), for 
instance for Delaunay triangulations (the graph then being its dual Voronoi diagram). 

This result further leads to a few very interesting questions: 


1. What is the behavior of the generating function of spin network evaluations when the edge couplings do not 
admit stationary points, and the stationary approximation fails? Normally we would expect an exponential 
behavior. Then, is it related to some exponential decay (e.g. of the 2-point function) in the non-critical Ising 
model in low or high temperature regimes? 


2 . As we pointed out above, the conditions (68) relating the dual triangulation to the couplings are scale- 
independent. They only depend on the angles and we can rescale all the edge lengths by an arbitrary factor 
without affecting the couplings Y e . Thus we actually have a whole line of stationary points by rescaling arbi¬ 
trarily the spins j e ’s. In |14j in the case of two-vertex graphs, it was found that this line of stationary points 
is related to the singularities of the generating function of spin network evaluations, i.e. to the set of zeroes 
of P r = Eccr lU c Y e . Can we expect also such a relation here? This requires a stationary phase evaluation 
of the generating function Z spin . Note that the zeroes of Pr are those of Z Ismg . Since Pr is polynomial (in 
Y e ) for all finite graphs T, the zeroes completely determine the partition function. In the context of statistical 
mechanics and in particular for the Ising model, they are called Fisher’s zeroes and their distribution in the 
thermodynamic limit determines the critical properties |39| . Could it then be that critical properties of the Ising 
model can be extracted from the stationary phase approximation of the spin network generating function? 


3. We could modify the generating function of spin network evaluations in order to get a single non-scalable 
stationary point, which would depend on the edge coupling Y e . This is actually the behavior of coherent spin 


below, these define a new generating function for spin network evaluations with a slightly different statistical 
weight depending on the spins j e . Although these other coherent states have a very nice geometric interpretation 
as semi-classical geometries, what would be their counter-part in terms of Ising model? 


network states introduced in [401 in the context of loop quantum gravity. As we will describe in section VIB 


4. Finally there is a clash with the correspondence with the Ising model. Indeed here the couplings Y e are allowed 
to run over all real (positive) values, while the correspondence with the Ising model requires bounded values for 
the Y’s given by Y e = tanh y e . This will be addressed in section EE below by generalizing the Ising model to 
O(n) models. 


B. Mapping Ising Correlations to Spin Averages 

The relationship between the partition function of the 2D Ising model and the generating function of spin networks 
induces a correspondence between the observables of both models. Those observables are on the one hand the Ising 
spin correlation functions, 


, T; i er v 2 


’ °w) = 


1 


Rising 


E 


®V2 ’ * ’ ^ 


E e y e v s (e)crt(e ) 


and on the other hand the expectations of products of colors, 

1 

in 

col 


(os? • ■ ■#> = • • -ft «(r, be}) n( tanh ^) 2j '- 


(70) 


(71) 
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where 2 j e £ N is the color of the edge e. The correspondence is obtained by taking derivatives of the fundamental 
equality, 


r^Spin ^glsing^2 n cosh 2 y e , \nZ Spin + 2\nZ Ising = ^21ncosh?/ e + 2(#V)ln2. (72) 

e e 

The two subtleties are the non-zero term on the right hand side, but it will disappear upon differentiating twice with 
respect two different edge couplings y e , and the dependence of the generating function Z Spm on tanhy e instead of 
simply y e like the Ising partition function, but this is easily accounted for by: 

Y = tanh y => dyf = cosh 2 yd y f . 


Taking a first derivative with respect to a given edge coupling y e gives: 

-- 77 —9 tan h Ve In Z spin = tanh y e - d Vc In Z Isin Y 

2 cosh y e 


(73) 


This is easily interpreted by introducing the nearest-neighbor correlations, i.e. correlations between the two Ising 
spins incident to the edge e, 


9 e = (o- s (e)CT t(e) ) = d ye In Z Ism9 , 


and the expectation of the color on the edge e, 

(2 j e ) = tanh y e dtanhy, In Z spm . 


Then equation (73) simply relates g e to (j e ), 

(j e ) = sinh y e (sinh y e - cosh y e g e ), 


s{e)&t(e)) — tanh y e 


1 


sinh 2 y ( 


(2 je) ■ 


(74) 

(75) 


(76) 


providing an expression for the mean color of an edge at fixed couplings, in terms of the spin-spin correlation of this 
edge in the Ising model. 

We then differentiate successively with respect to different edge couplings y e along a path of edges between two 
vertices in order to obtain longer range Ising correlations. For instance, let us start with three vertices Do, 1,2 linked 
successively by the edges e\ and e 2 . We differentiate with respect to y\ and y 2 and get: 


dyA f (HZ Isin9 )) 


dy e dy f Z ISinP 

Rising 


dy B Z ISin Bdy f Z ISin 3 

(Zi sin g )2 


(oo< 7 2 ) - (owl)( cti cr 2 ) 


(77) 


Y e Y f d Vc d Yf (In (Z spm )) = 


Y e Y f d Yt . d Yf Z spin Y e d Yc Z spin Y f dy f Z spin 


2 ^spin 


^£spin ^2 


= <(2 j 1 )(2 j 2 )> - (2 Jl )(2 j 2 ) 


(78) 


M - = si , lh2gi ( inh2i)2 (« 2 ji )( 2 h)) - 2 >) ■ ( 79 ) 

More generally, we consider two vertices on the graph T, that is an initial vertex vq = v and a final vertex v n = w 
linked by a path V consisting in n edges, e\ to e n . We consider cuts of this path: the cut P of the path V is defined 
by p intermediate vertices on V numbered v l P , i = 1 ..(p — 1), and divides the path V into (p + 1) smaller paths 
enumerated as Vp with I = l..(p+ 1). Differentiating the logarithm of the partition function and generating function 
with respect to the n variables yi to y n yields a sum over all such cuts P of the path V: 


WvVw)^ = 


—2 ra-1 

Ileev sinh (2je) 


(I[( 2 j e)) p, 

e£V 


(80) 


where we define the “connected” correlations as: 

p+i 

(a v a w )^ = ^(-l)^^)..^^), <I]( 2 je))?' = E(-!) P II< II ( 2 ^)) ■ 

P\V eer P\V 1=1 eG P : p 


( 81 ) 
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C. Distribution of the Edge Color 


Now we restrict our attention to the observables on a single edge, namely ((2 j e ) n ), for which the correspondence 
with the Ising model leads to explicit expressions. 

Theorem V.l. The exponential generating function of the spin averages ((2 j e ) n ) is 


Je(t ) = £«2? e n 
n> 0 


r 

n! 


1 

(l-0'e)(e‘-l)) 2 ' 


(82) 


It can be interpreted as the moment generating function j e {t) = X)n>o -P(2 j e = n ) e nt f or the following distribution, 


P{2j e = n ) 


n+ 1 ( {je) \ 

(l + 0e)) 2 \l + (je)J 


(83) 


where ( j e ) is given in (76). 


Before proceeding to the proof, we emphasize that this is not strictly speaking a probability distribution, since ( j e ) 
can be negative as we will see in the next section. This can be traced back to the fact that we are considering spin 
averages weighted by the spin network evaluations (which can be negative), which is not the expectation values of 
spin operators on coherent states, as underlined earlier in section [V A| 


Proof. Before beginning the proof per se, we want to emphasize that the simplicity of the result relies on fact that 
the derivatives of the Ising free energy In Z Isin9 with respect to a given coupling y e are all simple functions of the 
nearest-neighbor correlation. This is due to the fact that = 1 for all Ising spins. In particular, this implies 
Q^Z Ism9 /Z Ismg = 1, so that the second derivative of the free energy, which also is the derivative of the nearest- 
neighbor correlation, reads 


, , . dl Z Ising 

= dl In Z Ism9 = PP-r-. - 

ye Vising 


dy e Z ISin9 

Rising 



(84) 


The first thing to do is to relate the expectation ((2 j e ) n ) to the derivatives of the spin network free energy In Z spm . 
We start with the standard expansion of (2 j) n in terms of falling factorials, 

n 

((2j e ) n ) =J2 s (n,k)((2j e )^), (85) 

fc=l 


where x— = x{x — 1) • ■ ■ (x — n + 1) is the falling factorial and S(n, k) are the Stirling numbers of the second kind. 
Notice that 


<(2je)*> 


tanh y e y S pin 

gspin a tanh y c ^ 


( 86 ) 


The derivatives of Z spin can be related to those of In Z spm using Faa Di Bruno’s formula^] but it is simpler to use 
generating functions. Defining the shorthand notations 


it is well-known that if 


Tu = (d tanhVc Z spm )/Z spm and K n = 0? anh 


_^ £n 

m{t) = y 

' n\ 

n> 0 


then n{t) = 


E 

n> 1 


^n-r = 
n\ 


In Z spm , 

(88) 

In m(t). 

(89) 


The strategy is therefore to first find K n using the correspondence with the Ising model, then take the exponential of 
k (t) to get m{t) and extract its coefficients /i n . We prove by induction that 


K n = 2(n~- 1)! («i/2)” . 


(90) 


1 This gives here 


OV; In Z spin = 
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E (-i) 1+E " p 
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Ep PVp=n 
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This clearly holds true for n = 1. In addition, we recall from the definition of n n and the equations (75) and 

Ki = (2j e )/ tanh y e = 2 cosh 2 y e (tanhy e - g e ). 

Then the following holds: 

fi'n+l = ^tanh y e ^n = .! (ftl/2) ^tanhy e (^l/2), 

where the last equality makes use of the induction hypothesis. We therefore have to evaluate 

^2 = ^tanhy e ^l = 2<9tanh y e (cosh. 2/ e (tanhi/ e 9e)) 

= 4cosli' ! y e sinh y e (tanh y e - g e ) + 2cosh 2 y e ( 1 - cosh 2 y e d Vc g e ) 

= 2 cosh 4 y e (tanhy e - g e ) 2 = 2 (ki/2) 2 . 


(91) 

(92) 

(93) 


(using (9tanh y e = cosh 2 y e d Vl ,). From the second line to the third, we have used (84) and some simple algebraic 


manipulations. Plugging this expression into ( |92[ ) proves the formula for n n . The generating function n(t) is then 

(94) 


n> 1 




so that 


i(t) = 


— P K (*) 


(i-¥) 


2 ’ 


and fx n = [; t n /n\]m(t ) = (n + 1)! (^i/2) n . 


(95) 


Up to a factor tanh" y e , this is the expectation of (2 j e )—. We can now form the generating function of the expectations 
((2je)"), 


n>0 


i(*) = E<( 2 ^)") n! = 1 + EE S(n , /c)(/c + l)!(tanhy e K 1 /2) fe —- 

n>l/c=l 

j.n 

= ! + E(E S ( n ’ fc ))d)( fc + l)!(taiili2/ e Ati/2) 


(96) 


fe>l n>fe 


One recognizes the exponential generating function of the Stirling numbers of the second kind, J2n>k S(n, k)t n /n\ = 
(e* — l) fc /fc!. Therefore, together with (2 j e ) = tanhj/ e Ki, 

1 


j(f) = E(fc + 1 )(Oe)(e t -l)) = 


k> 0 


(l-0e)(e*-l )Y 


(97) 


In order to identify a discrete probability distribution, this expression has to be expanded onto powers of e 4 , which is 
readily done, 


«‘ ) = (rrWS ( ” +1) 


(je) 

1 + (je) 


and the coefficient of e nt is interpreted as the probability P(2j e = n). 

We can furthermore write the expectations as polynomials of order n in the mean color (j e ), 

n 

((2je)") =^2 s (n,k) (k + l)\(j e )\ 


(98) 

VT 


(99) 


fe =i 


or as 


<(2 jeD = E P ( 2 je = fc )^ 


fc >0 


(l + (je)) 2 


Li— n —1 


Oe) 


1 + (je) 


Li_-r 


(je) 


1 + (je) 


where Li„( 2 ) = J2k >l ^ * s the polylogarithm. 


(100) 
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D. 0(n) Models, Critical Ising Model and Phase Diagram 

Let us look at the potential critical behavior of the spin network generating function from the point of view of 
its duality with the 2D Ising model. For the sake of simplicity, we restrict our attention to homogeneous couplings, 
Y e = Y, Ve £ E. While exploring the range of all edge couplings is desirable from the point of view of spin networks 
and quantum geometry states, we see that it is not possible within the frame of the Ising model since Y is restricted 
to be smaller than 1, Y = tanhy < 1. To explore the regime Y > 1, we match the spin network generating function 
with the squared inverse partition function of the 0 ( 1 ) model, 

Z^n{Y) = l/(Z OW (Y)) 2 . ( 101 ) 

The 0{n) model, with n an integer, is defined as follows. Some n-component spins S) sit on the vertices of the graph 
and each equipped with a normalized measure dS such that J dS = 5^ v . The partition function reads 

Zo W (r) = / n dS V ll(l + YS sie yS t(e) ) (102) 

^ v e 

By expanding the product of the edge weights and performing the integral, a formulation as a sum over loop config¬ 
urations is obtained 


Z 0{n) {Y) = Y J ^ Y * E{1) ( 103 ) 

7 

where 0 ( 7 ) is the number of connected components of the loop configuration 7 and #E( 7 ) the total number of edges 
it covers. 

In the scaling limit, the 0{n) model gives rise for — 2 < n < 2 to the following the phase diagram (adapted from 

mi 

Y = 0 Y c Y 0 y = 00 

•--•-►-•--• 


Introducing 7 £ [0, 7 r) such that n = 2 cos 777 , it is found that pl^lflT| 

- Y = Y c is critical, with central charge c = 1 — 6 ^y^-, for g = 1 + 7 . Moreover Y c = 1/\J2 + \J2 — n on the 
hexagonal lattice. This is a second order phase transition with dilute loops. 

The region above Y c , i.e. Y > Y c , is called the dense phase and flows towards Y 0 whose location is Y 0 = 
1/^/2-y/2 — n on the hexagonal lattice. 

At Y = 00 , one gets the fully packed 0(n) loop model (whose universality class is lattice dependent). 

Our interest is the case n = 1, 7 = 1/3. Then, Y = 0 is the infinite temperature Ising model, with no loops at all. 
Y = Y c = 1 /a/ 3 is the Ising critical point with central charge c = 1/2, while lo = Y corresponds to the Ising model 
at zero temperature. 

We recall the expression which was found for the average spin of an edge in the Ising case, i.e. Y = tanhy, 

(je) = sinh 2 y — g e (y) sinh y cosh y, 


where g e {y) is the nearest-neighbor correlation at coupling y. Notice that the expectation of j e + 1/2 takes the simple 
form 


& + h = 

At small coupling Y = tanh y ■ 
the expected spin behaves like 


\ [cosh 2y-g e {y) sinh 2y\ e 2y (l - g e (y)) 


sinhy, g e (y) goes to zero as a power law (at 


+ e 2y (l + g e (y)) 


(104) 


least like y 2 if there is no 2 -cycle), and 


{j e )~Y{Y-g e (Y)), 


( 105 ) 


which obviously goes to zero, and thus (j e + 1 / 2 ) goes to 1 / 2 . 
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FIG. 5. The two plots on the left are plots of (j + 1/2} as a function of y on the hexagonal lattice. The leftmost plot shows 
the behavior from y = 0 to 1.7 (it goes to 0 at infinity). The phase transition takes place at y c ~ .658 around which the middle 
plot is centered. The rightmost plot shows the logarithmic singularity of the derivative of (j) at y c . 


When Y goes to Y 0 = 1, i.e. y — > oo, the correlation g e (y) goes to 1, and therefore the contribution e~ 2y (l + gJy )) 
goes to 2e~ 2y . To evaluate the behavior of 1 — g e (y ), one uses the low temperature expansion of the Ising nroder] it 
is easily checked that 1 — g e {y) decays at least like e~ iy . Therefore, 

(106) 

which also goes to zero. 

Furthermore, explicit formula exist for the the nearest-neighbor correlations of the Ising model. For instance, on 
the hexagonal lattice with isotropic coupling Y. one has (suppressing the edge dependence) [T?, '35] 

a(y) = coth(2 L(y)) (, a(k(y))A(y ) - b(k(y))B(y )), (107) 


with L(y) = \ In and k{y) = l/(sinh2 L(y) sinh2y), and 

a(k) = - ((1 + k)E(k i) + (1 - k)K(ki)) (108) 

b(k) = -(l-k)K(k 1 ). (109) 

7r 

Here K(k),E(k) are respectively the complete elliptic integrals of the first and second kinds of modulus k, and 
k\ = 2y/k/(\ + k ). Finally, 

A(y) = F 1 (arctansinh 2 L( 2 /), 1 — k 2 (y)) (HO) 

B(y) = ——— 1 ——,- ^F(arctansinh2T(y), 1 — k 2 {y)) — £ l (arctansinh2T(t/), 1 — k 2 (y))^ , ( 111 ) 


where k), E(<j>, k ) are the incomplete elliptic integrals of the first and second kinds of modulus k. 

The quantity (2 j) as well as its derivative with respect to the temperature can be rewritten in terms of standard 
thermodynamical quantities. Let us rescale y by the inverse temperature f) explicitly. Then the internal energy U ( (3y ) 
is such that 


Further, the derivative of (2 j) reads 


g(Py) 


u(Py) 

y #E ' 


( 112 ) 


9 / 3 ( 2 j) = 2 y( sinh(2/3y) - cosh(2/3y) 


dp In Z Isiny \ 

V# E J 


— sinh (2 fly) 


dj In Z Isin s 

y# E 


(113) 


2 It is an expansion with respect to the two configurations where all Ising spins are aligned, Z Is%ri9 = 2e y ^(l + C 2 ^~ 427 + 0(e~ Qy )). 
Corrections like e~ 2py come from flipping some Ising spins such that p edges have opposite spins at their ends. There is no p = 0 
correction since that would require a single edge to have opposite spins on its vertices and this is impossible if F is bridgeless. Then, 
C 2 is the number of pairs of 2-cut edges (i.e. such that cutting 2 edges disconnects T). Similarly, (7 i (J 2 ^~ Je yeCT s(e) cr t(e) — 

2e^# E (l + C' 2 e~ Ay + 0(e~ 6y )) with C' 2 < C 2 . 
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and introducing the heat capacity C{f3y) = /3 2 d| In Z Ismg , we get 

9p{2j) = 2y (smh(2j3y) + cosh(2/3 y) ") - sinh(2( 114 ) 

V v # E / p v # E 

Plots of (j + 1/2) are given in the figure [5j using the above formula for g(y). The Ising phase transition occurs at 
tanhj/ c = l/-\/3 (hence y c ~ .658) which corresponds to k(y c ) = 1. The singularity comes from the k ~ 1 behavior of 
a(k),b(k), as for instance 


b{k) 


k~ l 


1 -fc 2 

7T 


In 


16 

1 — fe 2 | ' 


In particular 0H> one finds around t = 1 — T/T c close to zero 


(2 j) ~ Cst x t In \t\ + analytic, 


(115) 


which gives a singularity of the type In \t\ for dp(2j). 

To conclude this section, the duality of the spin network evaluations with the 2D Ising model allows to identify a 
phase transition in the behavior of their generating function and corresponding coherent spin network state. This is 
clearly visible in the behavior of the average of the edge spin (j e ), which we computed from the 2-point correlation 
of the Ising model between nearest neighbors, which shows a discontinuity in its derivative. Let us point out that 
this average can be negative: as explained earlier in section V A here (j e ) is not the expectation value of the spin 
operator on the coherent spin network state but the average of the spin j e weighted by the spin network evaluations 
(which can be negative and do not define a true probability distribution). In this context, ( j e ) can not actually be 
interpreted as the average edge length of a dual triangulation to our graph, as suggested earlier in section [VA] by 
the analysis of the stationary point approximation of the expectation value ( j e )coh■ However, it does plays the same 
role as the nearest-neighbors correlations in the Ising model. It would nevertheless be enlightening to understand the 
geometrical interpretation of (j e ), in particular from the perspective of this stationary point approximation of the 
spin network generating function and its geometrical interpretation; doing this could shed light on a potential relation 
to a phase transition and continuum limit of discrete geometry in quantum gravity. 

To this purpose, we conclude this paper by the introduction of slightly different coherent spin network states, and 
thus of generating function for spin network evaluations, as defined in the context of quantum gravity m , which 
admit a different (and maybe better behaved) stationary point with a clear geometric interpretation, and which could 
be relevant for future investigation of this phase transition. 


VI. VERTEX INTEGRALS AND COHERENT SPIN NETWORK STATES 

We introduce a reformulation of the Gaussian integral for the spin network generating function over half-edge 
complex variables as a Gaussian integral over complex variables living at the vertices of the graph T. This will allow 
us to define new coherent spin network states with an improved behavior for the stationary point analysis. But it 
should also allow us in future investigations to tackle the case of graphs with nodes of arbitrary valency. 


A. The Generating Function as a Vertex Integral 


We use the spinor notations introduced in [36] and defined earlier in section V A We recall the e xpression of the 
generating function for spin network evaluations in terms of the angle couplings X a (Theorem III.5): 


z s P m ({x a }) = [ TT 

j C 4 # E tt 


„ —(z”| z”) , , ... 

e d e £e« (e) b) (e) ] g£ Q X a [z BM \z^ a) ) 


(116) 


where we changed the sign in front of y (z|^ by a linear symmetry sending ~z v e — > —z v z and having a trivial 

Jacobian equal to 1. This formula suggests considering angle couplings of the type X a = [Cs(a)ICt(a)) where Q! are 
arbitrary fixed spinors living on the half-edges similarly to the integration variables z v e . In that case, we can apply 
the following proposition, developed in the spinor formulation of loop quantum gravity: 

Proposition VI.1. Let Ll = (1, 0), Q,Zi,i € I be vectors in C 2 and I be a finite set; then the following holds: 



e -Gb) d 4 z 


7T 
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e Ei e /[Ci|2)<f2|zi)+[Ci|z][n|zi) 


elWfilGINIv). 


(117) 
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Let us point out that due to the anti-symmetry of the scalar products [C* 10 ) and l z i I z :i)■> the condition i < j is not 
particularly relevant. It only serves to avoid over-counting and does not reflect the necessity of a linear ordering. This 
proposition can be proved by putting together the following beautiful integral over SU(2), 

Lemma VI.2. 


dg&M’M = EjjpVl)! < 118 ) 

with the following relation between the Haar measure on SU(2) and the Gaussian measure over spinors: 

Lemma VI.3. From \ lf\ f 3&j . the integral of a homogeneous polynomial P{g ) in g £ SU(2) of even degree 2 J can be 
expressed as a Gaussian integral over C 2 : 




dgP{g) = 


R -(z\z) d A z 


(J + 1)! J c 2 7r 2 


P(z), 


(119) 


with g = |z)(fl| + |;z][fi|, with unit spinor fl = (1 0). 

The above results allow to reformulate the spin network generating function, introducing intermediate spinor vari¬ 
ables f v £ C 2 , living at the vertices v: 


Z S *' i "(r,{^a = [C. ( a)IC t (a))}) = 


/ C 4 # E J C 2 # v 


TT e {zJzJ d 4 z v e tt e e i;„ ie3 nciU W>HCIG][n|<> 

X X 77-2 X X 77-2 


We now use the following beautifully simple lemma of complex integration: 


( 120 ) 


Lemma VI.4. Considering an arbitrary holomorphic function f of spinors z £ C 2 , the following integration identity 
using the Gaussian measure over the space of spinors holds: 



e -(z\z) d A z 


f( z )e 


<*!€> 


/( 0 - 


This allows to integrate over all the half-edge spinors zf , leaving the generating function as an integral over vertex 
spinors only: 


Z S *> in (T,{X a = [Cs(a) |Ci(a))}) 



e (GIG) d 4 f v 


e -Ee[C (e, |(lG(e)><J t (c)l + IG(c)]K t (e)l)lC7 e> > _ 


( 121 ) 


This is the spinfoam amplitude, recently considered in [14j [lQ. These few steps allow us to bridge between the 
formulas for the spin network generating function as an integral over variables living on vertices as in mm and on 
half-edges as in m- 


This naturally leads to two questions. First, at a technical level, this vertex integral reformulation assumes the 
definition of the angle couplings from half-edge spinors as X a = [Cs(a) ICt(a))- This is different from the couplings 
defined from the Ising edge couplings, X a = -\/Y s ^Y t ^ a y How compatible are those two definitions? Since we 
already know that edge c oupli ngs are equivalent to angle couplings from the perspective of loop observables as 
explained earlier in Section IIC we can focus on asking how generic is the form X a = [Cs(a) ICt(a)) for angle couplings. 
We claim that it is completely general. Indeed around a 3-valent vertex, with angle couplings X\ 2 , X 23 , X 31 in C, 
there always exist three spinors |£i), IC 2 ), IC 3 ) £ (C 2 ) 3 such that X iti+ i = [CilCi+i)- Indeed possible choices are: 


ICi) = 




IC3) 


X 23 /X 12 

-X31 


Then we can act on these three spinors with an arbitrary SL(2, C) matrix without changing their products [Ci, Cj)- 
The second question is more intricate. This expression of the spin network generating function as an integral over 
vertex variables hints towards a similar reformulation of the Ising model possibly as an integral over odd-Grassmann 
variables living on vertices. We can also wonder if there exists a supersymmetry relating the new vertex variables f v 
to the Ising integration variables We leave these very interesting issues for future investigation. 
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B. Another Class of Coherent States and Generating Function 


We have reformulated the spin network generating function alternatively as a complex Gaussian integral over spinors 
associated to half-edges of the graph (Theorem |III. 5 ) or associated to vertices (Formula (1211), for angle couplings 
defined as X a = [C s (a)ICt(a)) i n terms of fixed half-edge spinors Q: 


z s P m (r, {x a }) = [ n 

^C 4#E %v 

1 n 


e <^l d 4 z 1 ’ j 2 e ( z sw\ Z ‘W] p J2 a x a [z s(oi) \z t(a) > 


e < g "l ^ d% j: eK (e)| (|G(e))(Ct {e) l + IG ( e)]Kt(e)l) IC* (e) > 


/c 2 # v ' 


We can generalize this re-expression to the spin network state defined as a gauge-invariant function on SU(2)# E /SU(2)# V : 


</>r,{A- Q }(3e) = / j f 


J c 2 # v 
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zV) V 

e d Z e p E e « (e) lff=ke (e) l P £c, *«[*.(*) I*t(a)> 


e -J 2 e [C ie) \H„ M g e Ht (e) |f^>> 


where we have defined the 2x2 matrices H v 's in terms of the integration spinors £„’s as: 

h v = + |£„][fi|. 


( 122 ) 


This expression naturally suggests to switch to SU(2) integrations instead of Gaussian integrals over the spinor space 
C 2 and we introduce the other class of gauge-invariant coherent states: 


br,{x a } 


(ffe) = TT dh v 

J SU(2)# V v 
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These are actually the coherent spin network states for loop quantum gravity and spinfoam models, introduced earlier 
in |40j , further developed in |14] and recently used in studying the coarse-graining of the spinfoam path integral in 
[15] . As one can see, the difference with the coherent states defined in Equation ( 60 ) and the above is that we have 
been using is the extra-factors (J + 1)!” 1 , turning the exponentials into Bessel functions. The two coherent states are 
actually related to each other by a Borel transform or reversely by a inverse Laplace transform. 

These new coherent states were introduced in the context of loop quantum gravity (in 3+1 dimensions) as describing 
good semi-classical states of 3d geometry. Here the geometrical meaning of those states is different since we seek for an 
interpretation in terms of 2D triangulations dual to our planar graph T. We will show below that they define slightly 
different weights for the generating function of spin network evaluations, which leads to a non-scalable stationary 
point uniquely fixed by the angle couplings X a , or equivalently the fixed half-edge spinor data Q. The introduction 
of the (J + l)! -1 factor for every vertex leads to a modified generating functional, to be compared with (351: 
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and the associated modified probability distribution for the edge spins, to be compared with (65): 
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(125) 


Here the edge couplings Y e are given in terms of the angle couplings X a according to Equation (19): 
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where e\ and e-i are the other two edges attached to the source vertex s(e), while ei and e 2 are the other two edges 
attached to the target vertex t(e). If we were to plug back the relation X a = ^/Vs(a)^t(a) i n that definition of the F’s, 
then it would be simply a consistency check. But here the X’s are defined instead in terms of the half-edge spinor 
variables Q. 

The stationary point analysis can be done exactly as earlier in Section |V A| Assuming that the spins f s are large 
and using the Stirling approximation for the factorials, we get the conditions at leading order in the spins: 


^14^ Js(e)(Js(e) ^jei){Js(e) 2j e2 ) Jt(e){^t(e) 2je 2 ) 

(Js(e) - 2 j e ) (Jt(e) ~ 2 je) 


(127) 


These equations are solved again if the spins j e are related to the edge lengths l e = 2 j e of a dual triangulation to our 
planar graph T with the condition that the edge couplings Y e are related to the geometric data by: 
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(128) 


in terms of the half-perimeters L v and the opposite angles 7 ". Unlike before, the Y e scale with the triangulation 
lengths and are not scale-invariant anymore. Given admissible couplings admitting a stationary point, we can not 
rescale the spins of that stationary point: we thus have an actual stationary point and not a stationary line as in 
Section IV Al 

Having started with the spinor da ta (Q 1 allows a finer analysis of the geometry of the fixed point. Indeed, comparing 
the two equations above (127) and (126), the conditions on the edge couplings for the stationary point are solved if 


for every angle a around every vertex v, we have the matching: 


J v {Jv 2 jot ) — |X a | — [Cs(a) ICt(a)) 


(129) 


where we have written a for the edge opposite to the angle a around the vertex v. More explicitly, around every 

(3-valent) vertex v, the spins j e at the fixed point are given in terms of the spinor data Q by the relations: 

J = ji+j 2 +j 3 , J(J-2ji) = ICCalCs)! 2 , ... (130) 

Following the work on coherent intertwiners in the context of loop quantum gravity and its discrete geometry [36l l47| , 
we can easily solve such equations to get the spins j e by introducing the 3-vectors V € M 3 corresponding to the spinors 
£ € C 2 by projecting them onto the Pauli matrices: 

v = \mo, \v\ = (cio, (131) 

where the Pauli matrices cr a , a = 1..3, are normalized such that they square to the identity. The previous conditions 
can be entirely re-written in terms of those 3-vectors: 

J(J-2j 1 ) = ^(\V 2 \\V 3 \-V 2 -V 3 ). 

Then, letting V) := |V)| and summing over all 3 such equations around the vertex, we define the total norm V = V) 
and the closure vector V = )Tb Vi, to get: 

4J 2 = (V 2 - V • V ), 2 ji = 1 _ _ (VV -V-V), (132) 

vV 2 -v-v 

where the spins at the stationary point are indeed fixed entirely by the value of the generating function couplings Q. 
If we further assume closure constraints on the initial spinor data, as introduced in [157} 140} , 

T; Vi = 0, or equivalently ^ ICi)(Ci| a H (133) 

i i 

the expression of the spins at the stationary point further simplifies to: 


2 J = V, 2 ji = Vi . 
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We obviously have to require norm-matching of the spinors on both ends of every edge e, in order to get the same 
spin j e from the equations at both it source and target vertices. 

To summarize the previous analysis, we start with spinor data living on the half-edges and satisfying closure 
constraints at the vertices and norm-matching on the edges: 

£lCXCI«I. Ve, (Cs(e) ICs(e)} = (Ct(e) ICi(e)) * (134) 

e3v 

In this twisted geometry setting [36] [371 i4S SS] , there exists a triangulation dual to our planar 3-valent graph such that 
the edge lengths U are the norm squared V z of the corresponding spinors. Then the new spin probability distribution 
p is peaked on a stationary point, such that the spins ji are half those edge lengths. 

This modified spin probability distribution has a very nice geometrical interpretation, which made it very useful 
to construct semi-classical states peaked on classical discrete geometries in the framework of loop quantum gravity. 
From our perspective here, the existence of the stationary point instead of the scale-invariant stationary line suggests 
that the corresponding spin network generating function may not have poles and therefore may not exhibit critical 
behavior. It would be very interesting to investigate this further and analyze in details the behavior of both spin 
network generating functions by a careful stationary point analysis involving not only the spin statistical weight but 
also the asymptotical behavior of the spin network evaluations (as studied for instance in e.g. m), and study the 
counterpart of those from the viewpoint of the 2D Ising model. 


VII. CONCLUSION & OUTLOOK 

In the present paper we related the partition function of the Ising model on a planar trivalent graph T to the gener¬ 
ating series of the spin network evaluations on the same graph. Although we explored some of the first consequences 
of this duality, many open questions remain which we will explore in the future. 

• Is it possible to extend the correspondence to non-planar graphs and/or graphs with vertices of arbitrary degrees? 
In |l2j a general formula was obtained for the trivalent non-planar case. In loop quantum gravity, coherent states 
on graphs with vertices of degrees greater than 3 also exist. We expect that those formulae to be related on the 
Ising side to formulae like the Kac-Ward formula (see !49]), where Kasteleyn orientations are extended to spin 
structures. 

• In jT3] a Gaussian integral formulation was provided also for spin network evaluations of graphs equipped with 
SL 2 (C) holonomies on the edges. From the point of view of Loop Quantum Gravity, they are natural objects 
to study. Is there a counterpart of holonomies on the side of the Ising model? What is the link with the spin 
network evaluations? 

• When the graph is embedded in the 3-sphere or 3-hyperboloid, spin network evaluations have famous deforma¬ 
tions known as quantum spin networks , based on the q-deformed ?? 9 (SU(2)). Is there a statistical model which 
maps to the quantum evaluations via supersymmetry and which deforms the Ising model? In the positive case, 
what is the interpretation in terms of this model of the conjectures on the asymptotic behavior of quantum spin 
networks (US:, US], [201)? 

• Our analysis in Section |V A| is a first analysis of the relations between stationary points of the distribution of 
the spins induced by spin network coherent states and critical values of the couplings for the Ising model. Can 
we better understand this relation, which at this stage might just look like a coincidence? Our analysis gives a 
relation with the critical temperatures for the Ising models only if the graph is isoradial. Can one push further 
the analysis to non-isoradial graphs? 

Then there is much to do on the interplay between spin network evaluations and the Ising model in the thermody¬ 
namic limit. In particular, we might be able to use methods developed for the Ising model to study the generating 
function of spin network evaluations. Also comparing the equations which appear in the resolution of both models 
(Biedenharn-Elliott identity and Yang-Baxter equations for instance) should lead to a fruitful cross-fertilization. This 
opens new perspectives on the coarse-graining of spin networks and their spin foam transition amplitudes for quan¬ 
tum gravity, which might meet the intertwiner dynamics models developed in m and their interplay with condensed 
matter models. Here is an outlook for potential future investigation: 

• The fact that the stationary points found in Section [V A| for the distribution induced by the coherent spin network 
state on the length operator are the same as the critical values for the Ising model in the thermodynamic limit 
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seems a bit mysterious. The stationary points we found should be compared with those of the stationary phase 
approximation of the generating function of spin network evaluations Z spm itself. As the latter takes the form 
Z spxn = 1/Pp where Pr is a polynomial, we expect the stationary phase method to reveal the singularities 
of Z spin . They obviously are the zeroes of Pr, hence the Fisher’s zeroes of the corresponding Ising partition 
function, as Z Ismg ~ Pp. This is a first way to study the thermodynamic limit, as the distribution of Fisher’s 
zeroes in this limit carries information on the critical behavior |39j . 

• Coarse-graining techniques are well-developed for the Ising model. The duality we exhibited suggests that 
they could be applied to spin network evaluations. We expect the star-triangle transformation to be directly 
the Ising equivalent of the behavior of spin network evaluations under Pachner moves. For instance, is there 
an equivalence between the Yang-Baxter equation for the Ising model and the Biedenharn-Elliott identity for 
Wigner symbols? 

• The duality of the 2D Ising model between low and high temperatures is expected to carry over to the spin 
network side. Does it then correspond to some UV-IR duality in 3D quantum gravity amplitudes? This typically 
requires to use non-trivalent graphs, as the duality pushes a model from the primal graph to its dual. 

• Finally, properties of spin network evaluations might be described in the continuum as a conformal field theory, 
dual to that of the continuum Ising model. Indeed, on the one hand, the critical 2D Ising model is conformally 
invariant and can be described as a Wess-Zumino-Witten model for the coset SU(2)! x SU(2)!/SU(2) 2 (were the 
indices denotes the level of the theory) (see e.g. [5T]1. One can consider its “inverse” by going to negative level 
or equivalently replacing the gauge group SU(2) by the coset SL(2,C)/SU(2) j55j. On the other hand, the spin 
network evaluation live on the boundary of the spinfoam model for 3D quantum gravity, which is described in 
the continuum by a Chern-Simons theory in the bulk for the Poincare group inducing a Wess-Zumino-Witten 
theory on its boundary. It would be enlightening if the two sides of the story would meet. Could this fit in the 
paradigm of gravity/CFT correspondence? 
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